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ABSTRACT 


A finite element method was applied to the unsteady 
transonic small disturbance equation and integrated until 
the solution converged to the steady state for a thin non- 
lifting airfoil. The method of weighted residuals was used 
to formulate the finite element equations, and Houbolt's 
method of central differencing in time was used to integrate 
these equations. 

A secondary investigation applied the steady transonic 


small disturbance equations to a converging-diverging nozzle. 
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EN TRODUCTTON 


Transonic inviscid flows past a smooth airfoil may be 
expressed in terms of the velocity potential ® satisfying 
the transonic small disturbance equation, 


(1 - M - МО + и, = 0 (1) 


This equation presents two major difficulties, 1) it is non- 
linear and 2) it is of the mixed hyperbolic-elliptic type. 
Analytical solutions to non-linear equations are difficult 
to obtain. One must normally resort to numerical methods. 
When the coefficient (1 - м“ - м“ 


со 


(ү + E) in equation 1 
is negative, the flow is supersonic and the equation is 
called hyperbolic; otherwise the flow is subsonic and the 
equation is elliptic. The forms of the two solutions are 
fundamentally different. Hyperbolic equations admit both 
О ыш ie which propagate only in characteristic 
een er ons, sand the presence of shock fronts. Elliptic 
equations, on the other hand, require that the dependent 
variables and their derivatives be continuous and that a 
ο πας ап part of the flow field affects every other 
part. Many non-linear elliptic equations are solved by 
appropriate relaxation iteration techniques by casting the 
equation in Poisson's form with the non-linearity actıng as 
Ai ОО ОК Осе solutions to hyperbolic equations are 


usually obtained by the method of characteristics or by 





A erence marching techniques which use an artificiali 
viscosity to represent the average jump conditions across the 
shock wave. In mixed supersonic and subsonic flows, normal 
numerical procedures break down because the boundary between 
the two regions is not known a priori. 

Finite element numerical techniques have evolved as a 
powerful tool in obtaining approximate solutions to a wide 
variety of engineering problems, particularly ones with 
Neumann-Dirichlet boundary conditions, i.e., elliptical prob- 
lems. They offer several outstanding advantages. Some of 
these are: 


1) Non-homogeneous problems may be treated with relative 
smolıicity. 


2) Complex geometries may be modeled with relative ease 
since the elements can be graded in size and shape to 
follow boundaries of arbitrary shape. 

3) Once the finite element model has been established, 

a variety of problems can be solved by supplying the 
computer with the appropriate data. 

Chan and Brashears [Ref. 5] developed a finite element 
computer program for steady transonic flow over a non-lifting 
airfoil. This program uses the least squares method of 
weighted residuals to approximate equation 1 by a system of 
algebraic equations, and assembles the equations in a special 
τπτ τος the hyperbolic region of flow. This tech- 
nique prevents the influence of downwind nodes from propagating 
Meee n the supersonic region. 

ο προς ο μες thesis is to investigate the possi- 


bility. of speeding the convergence of a solution by transform- 


Mme the steddy transonic equation to the unsteady equation 





cer ratine Over time Until the time dependent terms 
ἵ πια απο то extend the program of Ref. [5] to the tran- 


AO econ ош а Converging diverging nozzle. 


1n 





ΠΡ DISCUSSION ΠΕ iN IE ELEMENT APPROACH 


In a continuum problem of any dimension, the field 
variable, whether it is velocity potential, velocity, temper- 
ature, displacement or some other quantity, possesses infin- 
itely many values because it is a function of each generic 
pommes solution region. Consequently, the problem is 
oM οἳ απ intinite number of unknowns. The finite element 
approach subdivides the solution domain into a finite number 
of subdomains called elements and expresses the unknown field 
variable in terms of assumed approximating functions within 
each element. The approximating functions are sometimes 
called interpolating functions and are defined in terms of 
the values of the field variables at specified points called 
nodes or nodal points. Nodes usually lie on the element 
boundaries where adjacent elements are considered to be con- 
nected. In addition to boundary nodes, an element may also 
contain interior nodes. The nodal values of the field variable 
Paget iceimtorpOolating function for the elements completely 
define the bahavior of the field variable within the elements. 
Once these unknowns are found, the interpolating functions 
define the field variable throughout the assemblage of the 
element. 

clearly, sehe nature of the solution and the degree of 
accuracy of the approximation depends not only on the number 


and size of the elements used but also on the interpolating 
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functions selected.  Interpolating functions may not be 
chosen arbitrarily because certain compatibility conditions 
must be satisfied. Often such functions are selected so 
that the field variable or its derivatives are continuous 
across adjoining element boundaries. Once the problem is. 
formulated in terms of individual elements, the contributions 
of each element may be assembled to define the entire solu- 
tion domain. This means, bL example, that if we are treating 
a problem in stress analysis we can find the force-displace- 
ment or stiffness characteristics of each element and then 
assemble the elements to determine the stiffness of the whole 
structure. Finite element solutions are not, of course, re- 
stricted to structures problems, but the matrix of equations 
defined by the interpolating functions and the nodal field 
variables is still referred to as the stiffness matrix 
regardless of the field variable in the problem. 

ОШИО to Continuous problems by the finite element 
approach always follow a systematic step-by-step process. 
This process is completely general to the finite element 


method and it is outlined below. [Ref. 4] 


Еее the Continuum. 


The first step is to divide the solution domain into 
elements. A variety of element shapes may be used and one 
or more different element shapes may be used in the same 
region. The type and number of the elements used in a given 


problem are a matter of engineering judgement. 


7 








Be Sect thesimterpolatins functions. 


ὑπ πο Step 1s te Choose Che type of interpolating 
function to represent the variation of the field variable 
over the element. The field variable may be a scalar, a 
vector, or a higher order tensor. Often polynomials are 
selected as interpolating functions for they are easy to 
integrate and differentiate. The degree of polynomial 
chosen depends on the number of nodes and the nature and 
number of the unknowns and the continuity requirements 
imposed at the nodes and the element boundaries. The un- 
known quantities at the nodes may be assigned to the field 


variable and their derivatives. 


3. Find the element properties. 


After the elements and. their interpolating functions 
have been selected, the matrix of algebraic equations which 
express the properties of the individual elements must be 
determined. Several methods are available for this task. 
These are: the direct approach, the variational approach, 
the method of weighted residuals and the energy balance 
approach. Reference [4] is a good source of information 


on the various techniques. 


4. Assemble the element properties to obtain the system 


equations. 


ιν lesproperteies Of the over-all system, the matrix 


equations expressing the behavior of each element must be 


195 








added to the matrix equation of all other elements. The 
basis for this assembly procedure stems from the fact that 
connecting elements have common nodes and the field variable 
must be the same for each element sharing that node. 

At this point the boundary conditions for the system of 
equations must be accounted for and the system of equations 


must be modified before it is ready for solution. 


5. Solve the system of equations. 


The assembly process of step 4 produces a set of simul- 
taneous equations which can be solved to obtain the unknown 
field variables. Linear equations have a number of standard 
solution techniques readily available, but solutions to non- 


linear equations are more difficult to obtain. 


6. Make additional computations if desired. 


Important parameters, such as pressure coefficient in 
aerodynamic problems, may now be calculated from the values 


of the field variables. 
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III. THEORY AND BASIC EQUATIONS 


A. STEADY TRANSONIC FLOW 

Chan et al. [Ref. 5] developed an algorithm to analyze 
steady transonic flow over non-lifting thin airfoils.  Bound- 
ary layer effects were ignored and the imbedded shock wave 
was assumed to be weak. These assumptions are consistent 
πι παμςοπις small disturbance theory which can be expressed 


mathematically by the following expressions. 
(1 - μὲ - МАС + Dodo. м» = 0 
со co КОО Y (1) 


Boundary conditions - 


у. ф-= 0 at infinity (2) 

A de ας on the airfoil (5) 

u = 0 on the line of symmetry 
Le 


where g(x,y) defines the airfoil and dg/dx describes the 
ESO ST OD. 

. The above expressions appear in their dimensionless form 
where $ = perturbed velocity potential and the perturbed 


veloeityzceomponents in the x and Y directions are respectively 


defined as 


и = oe 
V = y 


15 





М, - Ireestream Mach number and y = ratio of specific heats 
ms που αιτ πο taken το Бе 1.4. The physical coordinates 
x' and y' and the velocity potential $' are related to the 


dimensionless quantities by 
Mes y. с $ = $'/cU. 
where c is the chordlength of the airfoil and Ii ene trees 


stream velocity. 


Once the flowfield solution is determined in terms of the 
perturbed velocities, the secondary unknowns are subsequently 


calculated. These include: 


U 1 
А 2 ο. 21 — 
а = [55002 - У) + (ql (5) 
Р 
O 29 Y/ (Y + 1) 
2 2 
E zu PNEU V "MU 
BUDE о | (8) 
со ο. ο. со 
Е ее — 11 5 normalized treestream velocity, 


a = local sound speed, p = local static pressure, M = local 
Mach number, V = total velocity, Po = stagnation pressure 


and С E s ως coefficient. 


B. UNSTEADY TRANSONIC FLOW 
Ui s l csnsaenuceenwiscrd Zlow may be expressed in terms 
of the velocity potential ¢(x,y,t) to a first order approxima- 
Cion Dy 
16 








τ (9) 


2 2 | 2 
"C RED sp) ch o EM ОУ cc 


уу 
which has the same non-linear coefficient retained for steady 
eransonTe Том in Equation 1. 
Boundary conditions require that the disturbances vanish 


БО гош the airfoil, 


$ = 0 Py πο» (10) 


and that the flow remain attached to the body. Let B(x,y,t)=0 
define the body at any instant. The surface tangency restraint 


may now be expressed by the substantial derivitive DB/DT 


vanishing. 


DB/DT - B, * (1 * 0,)B, * 6, B, (11) 


If the body remains stationary B. - 0, and the tangency 


condition becomes the same as in steady flow, namely 
v= (1 + u)dg/dx (12) 
where dg/dx represents the airfoil slope. 


The pressure coefficient for isentropic unsteady com- 


pressible flow is defined by 





с = — fu ασε μὲ (24 +24 + 92 + AD. 
UE á (13) 


Expanding by the binomial expansion and retaining only the 


first order terms gives, 


со - 2% (14) 


р t 


1% 





IV. FINITE ELEMENT FORMULATION 


The method of weighted residuals is a technique for 
approximating solutions to linear or non-linear partial dif- 
ferential equations and it is the basis for the finite element 
formulation of the transonic small disturbance equation 
(Equation 1). This procedure involves assuming the general 
functional behavior of the field variable which would approx- 
imately satisfy the basic equation and boundary conditions. 
substituting this approximation into the original differential 
equation results in some error called a residual. A system 
of algebraic equations results when a weighted average of the 
residual is forced to vanish as it is averaged over the entire 


domain. 


ΛΑΟ FLOW 


The approximate solution to equation l is assumed to be 


^ 


ф = М.ф. (15) 


where .N; are the interpolating functions which exhibits the 
behavior of equation 1 and $, are кше кипаегетттпеа рагатегет< 
at the nodal points. 

When s substituted into equation 1, the resulting 
residual is 


В = [1 - Μ΄ - MECH FLING IN xx + №. уу (16) 
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The weighted average is determined by multiplying the 
residual R by m linearly independent weighting functions W. 
and integrating over the elemental domain. Forcing this 


residual to vanish yields, 
SIW RdA = 0 = | Comm 


Chan et al. [Ref. 5] found that when the weighting function 
W. for equation 1 is chosen to be 9 R/ 99; τυ πι. 
is positive definite and well conditioned. This choice of 
weishting functions is referred to as the method of least 
squares because it is equivalent to minimizing the square of 
the residuals summed over the domain with respect to the un- 


determined parameters. That is, 


X f J R^dA 


i 


эх/э $ = ЈЈӘК/Ә Ф. КАА = 0 (17) 


Integrating over the domain produces the system of 


algebraic equations 


ην. Ξ. 0 18 
ij?j (18) 
where the elemental matrix 513 is defined as 


Si; " IF; dA (19) 


With 9; апа Р. equal to 


κ,χϑκ' Ni xx у № yy 


Q. = [1 - M“ - Mé (y + 1)N 


J 


Ὃ 
H 


2 
m E ud PO x 
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B. UNSTEADY FLOW 

Development of the unsteady flow finite element equations 
is similar to the procedure used to formulate the finite ele- 
ment equations for steady transonic flow. The least squares 
method of weighted residuals is again used but $ is now a 
function of time as well as the spatial coordinates x and y. 


The approximate solution has the form, 


ç = N, (x,y)#; (t) (20) 


Substituting ¿ in the unsteady transonic small disturbance 


equation, the weighted residual becomes, 


ES RO Rz) dA (21) 


where 


I 
li 


2 2 
ατα τ ΜΕ ΜΕ + LON ING xx * NG уу; 


2 9 
о uat; 


DÀ 
1 


Aj 
I 


MNS 
Expanding equation 21 yields 

x = SSIR + R? + RZ + 2R)R, * 2R,R;]dA (22) 
and minimizing with respect to the undetermined parameters 
b: the following system of algebraic equations result, 


$; 7 ' 2/ [9R 1/894 [R4 + К, + R¿]dA 
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where Р. has been previously defined in the steady finite 


element formulation 


The above equation may be rewritten in the form 
$..6. + SC.:d. + SM..6, = 0 0 23 
UE UE ij*; ο 


The stiffness matrix Si; is the same as that developed 
for the steady transonic equation and the damping ee and 


mass en matrices are defined below. 


SC.. = -SIMEN. _P.dA (24) 

nou uS x i 

SM.. = -SIMEN.P.dA 25, 
η Е 
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V. ELEMENT DESCRIPTION AND ASSEMBLY OF EQUATIONS 


ООО ELEMENT DESCRIPTION 

The basic element used in the finite element program is 
the non-conforming cubic triangular element developed by 
Bazeley et al. [Ref. 2]. Also used in the program are the 
ο ο πο σα] апа trapezoidal elements constructed from these 
triangular elements. These three types of elements can be 
mixed and used freely in the entire flow region except that 
only trapezoids should be used to cover the supersonic and 
mixed region in order to enact the special assembly procedures 
required by the hyperbolic equation which describes the flow 
MENA e region. 

The basic triangular element is shown in Fig. 1, which at 
each vertex has the velocity potential and the velocity com- 
ponents as the undetermined parameters. This type of element 
was chosen because both Dirichlet and Neumann boundary condi- 
tions can be imposed with equal convenience. In addition, 
because velocity components are used as primary unknowns 
secondary parameters, such as Mach number and pressure coef- 
ficient can be calculated directly without resorting to 
numerical differentiation, which would produce additional 
errors: 


InsrhezeTement, the approximate solution is assumed as 


ф = м, $. (i = 1 to 9) 
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3(9, us. Ms) 





Bene c leranguiar Elemente 


2:9 





In which q; 's are the nine undetermined parameters of $ and 
n are the interpolation functions which are expressed in 
terms of the area coordinates. 

The shape or interpolating functions and their first and 


second derivitives are defined below. 


ВЕТ, 
TEE IE SE 
Dg = a. 
НЕ 565 
Шс те От апп е l= 2-5 = (ο οκ - b4C,)/2 
cu CUIR S) (ск š ce) 
ου Ы; го) 
Н = 045 55% 


Hy = 6:5:5к * 6:65: * 6645) 


ы у= ү us Кот т C1045 


y be J 


Hex ^ 2(zkb;jbk + cjbkbk + cxkb;b;) 


Hyy = C к + СКС + Bue oq) 


with i 20008 7. 5) κ Ξ qo τηε. οποια’ for 1. = (1.4.7), 
ОЭ) 

m _ 
М, = z (3 20) + 2H 


Z 
lI 


[6b.c, (1 - 55) + 2H_ 1/24 


= 
II 


ESL LPS 
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= 2 _ 2 
Ni xx = [6b. (1 25.) + 2H >] EAS) 
N era: ооду 
куу SS: 1 yy 
for 1 = (2,5,8), 1 = (1,2,3) 


ИШ. c7 R 
N. = š; С τ τ ΠΗ 


P 


_ | 2 
ШЕ = [2b; 5; (>; ©, ο + 25. + S [ZÈ 
DN - [2€; 51 (e 9; - €; бк) * л A 
= 2 ] 2 
et [ 2b: CAS с: бк) + ДЫ, (2А) 2. + πα / (24) 
N = wes (45. - c,o,) + OH] / (2A)? 
J y i j jk уу 


for 1 = (3,6,9), 1 = (1,2,3) 


"m: . 
1 ` S£ (δ1ζκ - Ὀκξι) + BH 


= 
| 


= 
u 


[2b C; (b ty - 6х5.) + 8н] /24 


Rx 
= ü 2 
Ny y = 120,5; (6:5; - 65.) + 25 + 8H] /24 
N = re τν 5ης ΡΤ] ΠΠ. 
и хх 1 j K К”) хх 
200 | , 
№ yy EE Б bk.) * 4c1(2A); + PH A, 
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Quadrilateral and trapezoidal elements as shown in Fig. 2 
are also used in the present program, the former is used in 
the subsonic region and the latter in the mixed and super- 
sonic region. The element matrix for the quadrilateral ele- 
ment is obtained by combining appropriately the matrices for 
two triangles, while the matrix for trapezoidal element is 
obtained by averaging contributions from two left-running and 
two right-running triangles. The averaging process is designed 


to remove the bias effects inherent in the quadrilaterals used. 


B. ASSEMBLY OF EQUATIONS 

Straightforward application of the finite element assembly 
technique to transonic flows would fail (the solution diverges) 
Beeause nis would allow disturbances to propagate upwind in 
the supersonic region of flow where the governing equation is 
hyperbolic. Hyperbolic equations have a time-like dependency 
in that the solution at the downwind station is affected by 
the upwind station but not vice-versa. Assembly techniques 
for a transonic flow finite element program must take into 
account this time-like dependency. If the x-axis is taken 
as a time-like direction in the supersonic region, the ele- 
ment matrix may be assembled in a way similar to a backward 
finite difference operator, which has been successful in 
solving hyperbolic equations. 

Consider the rectangular element sketched below with the 
upwind station I and the downwind station Il, each having two 
ЕС рона with the element type chosen. The element 


matrices can be constructed in the usual manner. ` 
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а. 


Quadrilateral 
Element 


b. Trapezoidal 
F Ten 


Figure 2 - Quadrilateral and Trapezoidal Elements 
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However, before assembling the element matrix into the system 


Bere: ıx tnemen-linear coefficient of equation l is evaluated. 
С=1- М - м (ү + 1) и 


The sign of the coefficient being positive, zero, or 
Dee ıcive derines the equation as elliptic, parabolic, or 
hyperbolic. If C is non-positive for all nodes in the element, 
the rows representing the improper-downwind influence on the 
solution at an upwind station are ignored during assembly. 
This feature is automatically applied in the program requiring 
only a little care in arrangement of the nodes of the element. 
In the anticipated supersonic region, element node points 
should be arranged in the order as indicated in figure 5, 
starting with the upper left corner and proceeding in the 
Goumeer-clockwise direction. In the elliptic region, 1.e., 
where the coefficient is positive, no special assembly tech- 


nique is invoked. 


MEME D DIVE PROCEDURES 
With the equations assembled and the proper boundary con- 
ditions imposed, the system of non-linear algebraic equations 


is solved by iterative procedures in the form 


POE 
оо. E 5 


1 


to solve for the solution in the nth iteration. The function 
$ is defined as 


N N 
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FLOW 
DIRECTION 


ы с > 


Figure 3 - Nodal Arrangement for Supersonic Region 
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паса спе under-relaxation factor 9 is in the range 
EIL rorosubsonse flow 9 = 1, which is simply a succes- 
On approximation, yields good results, but it is necessary 
to under-relax somewhat with 9 approximately .5 for super- 
critical flow. Generally, a smaller relaxation factor will 
Meee solution more stable but it will tend to slow down 
πες ος convergence. 

Equation 23 is subject to the convergence criterion that 
the change in local Mach number between two successive iter- 
ations is less than a prescribed value € at all nodes in the 
Fiowserield: That is, 


. (n) a qM(n- 1) 
uei "ues 
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war INTEGRATION OF UNSTEADY FINITE ELEMENT EQUATIONS 


The unsteady transonic small disturbance equation (equation 
9) when suitably reduced to a finite element approximation 


appears in the form, 


- SO NON INE Ss 26 
ijj 1593 Mig? j j cay) 


This equation is analogous to a damped spring mass system, 


Mence S.., 
1] 


the stiffness, damping, and mass matrices. 


С апа = are respectively referred to as 


Mathematically, equation 26 represents a system of second 
order differential equations with constant coefficients, which 
can be solved by standard numerical procedures for differential 
equations, such as Runge-Kutta or Milne methods. However, 
this wuld be a very costly technique if the coefficient 
matrices are very large. In practical finite element analysis 
there are a few effective methods which take advantage of the 
banded matrices usually encountered in finite elements. One 
such method is the direct numerical integration method. 

Direct integration involves a numerical step-by-step 
procedure aimed at satisfying equation 26 only at discrete 
time intervals Оаа осие Dp тте с. Conceptually, 
direct integration is a finite element method in space and a 
finite difference method in time. Examples of direct inte- 
gration are the central difference method, Houbolt integration 


and the Wilson method. The first two schemes are finite 


Y 











difference schemes whereas the latter is a linear accelera- 
tion method. Linear acceleration integration assumes a 
Linear variation of acceleration from time t to time t -*At. 

Central differencing can be very effective in the solu- 
tion of many dynamic problems especially those that involve 
a large system of equations. However, this method is unstable 
for all time steps larger than a critical time step. 

Houbolt integration is an implicit finite differencing 
method related to central differencing, only it has the ad- 
Montage of being stable for all TIME STEPS. 

The Houbolt method was used to integrate the unsteady 
finite element equations because of this stability. 

Houbolt integration uses the following finite difference 


expansions: 


= И _ 2 
ΜΠ με ieat Si ь- 246] 1788 


[116 189 / ОА 


Pi trat t+At VE AT AN 


which are two backward-difference formulas with errors of 


order у 


” 


The solution of $i must satisfy equation 26 and at 


,tt+At 
time t+At equation 26 becomes 


21393 ,t«at* 701595 ceat * 935495 teat 


Substituting the finite difference formulas for E COME and 
3 


rearranging all known vectors on the right hand side, the 


Solution ror o. 


Р UNIS obtained, namely: 


52 





ф - 
EE uh ^u seat o By cent (27) 
MOL uua oe a SM. + aC ель 
Where the constant integration coefficients are: 
2 2 
am ЛЕ 
a, = ОЕ 
a, = Sort 
a, = τι 
a, = 28g 
a, = -a,/2 
ag 7 ag/ 2 
a, = a,/9 
Equation 27 may be written as 
ἘΞ 19; erat = КЕ, (28) 


where the effective stiffness matrix EE and the effective 


load vector = are defined as: 
/ 


FE q 


ij ij uL MD 


DI E EE °з? 


$ 


i j t-2At) 


ee 7 


Accurate knowledge of the vectors Е and ИЕ 


Be eregurned to yield an accurate solution for A COE and 
2 


EIS 





normally the Houbolt integration scheme requires a special 
ne procedure to determine the initial two vectors 

$5 At one ПЕ ceros рт тшату interest оғ 
this problem is to integrate the equations until they con- 
verge to a steady state solution, it is not necessary τυ 
obtain an accurate time history of the flow. Errors induced 
by the inaccurate starting vectors will vanish as time 
ΠΠ infinity. Therefore, the starting vectors may be 


chosen somewhat arbitrarily. 
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VII.  CONVERGING-DIVERGING NOZZLE 


Far Shen Ret'” 0] demonstrated the feasibility of 
calculating compressible flows through a converging-diverging 
(Laval) nozzle by dividing the region of calculations into 
ance patches, a subsonic region, a supersonic region and a 
transonic one, of course bounded by the other two regions. 
The locations of the boundaries for each region were chosen 
arbitrarily provided the sonic line is bracketed by the 
subsonic and supersonic boundaries. 

Two different finite element formulations were used for 
the subsonic and supersonic regions, but Shen [10] resorted 
to analytical approximations to cover the transonic patch. 
Ns restricted the calculations to nozzles with small throat 
curvatures because no analytic solutions exist for nozzles 
with large throat curvatures. It is conceivable that STRANL-II 
could be adapted to provide a continuous solution throughout 
all three regions. 

Outside the transonic region of flow the governing small 


eur panes equation is 


2 - 
T E ο ο ο (29) 


This holds for both subsonic and supersonic flow. Comparing 
equation 29 with equation 1, the transonic small disturbance 
Equation we notice that Only the non-linear coefficient 


Mé (y + l)u distinguishes the two equations from each other. 
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This coefficient becomes negligible when the Mach number 
weones less than 28 or greater than 1.2. With this con- 
sideration in mind, it was assumed that equation 1 would 
adequately describe the flow through the nozzle and that 
nte element formulation developed for the non-lifting 


airfoil would apply to the Laval nozzle. 


A. BOUNDARY CONDITIONS 

Two solutions are possible for a converging-diverging 
nozzle: 1) Symmetric flow, where the flow is subsonic through 
the domain, except for a small supersonic region near the wall 
in the throat, and 2) Asymmetric solution, where the flow 
accelerates to sonic velocity in the throat and then continues 
to accelerate to supersonic velocity in the diverging section. 
Different boundary conditions apply for the two solutions. 
Бере Symmetrie case, both inlet and exit velocities must 
be specified. Inlet and exit velocities are equivalent in 
the subsonic solution. The supersonic solution requires that 
only the inlet velocities be specified. If the exit veloci- 
ties are also applied, the problem is overspecified and the 

/ 

БОГЕТ ОП may not converge. 

Velocities at the inlet and exit are not uniform in the 
y direction, therefore the disturbances cannot be set to zero 
as in the case of the non-lifting airfoil. Boundary velocities 
must be calculated by solving equation 29 analytically. 

Equation 29 is a linear equation which can be mapped to 


πα s equation, 





by letting y' =, / (1- M^) XM eS оп пы еаз11у 


solved for the case of the hyperbolic nozzle by transforming 
from cartesian coordinates to nic eoordımates This 
transformation simplifies the solution because the stream 
Mnes must be hyperbolas to follow the nozzle boundary and 
therefore follow the hyperbolic coordinate v= constant. 

e ртс coordımares y and v are chosen such that 
the curves u - constant are ellipses and the v curves are 
, τι.’ then the velocity potential which satisfies 


acc scequation tor a hyperbolic nozzle is simply 
ф = Au 

where A.is a constant of integration. The stream function is 
ф = Av 


The transformation w = u + iv = cosh * (2z/a) gives rise 


to the elliptic coordinates 


ЕТА СОБЕ = 1/2 а sinhu sin v 
2 A 

т = 1/2 a [coshu - sin ν] 

в = AV + а/2)^ + x? 

Boc (y - a/2)* E x? 


Solving for y and v produces 


1 


Е το.” Mr a 


v= cos”! prora] 


The nozzle boundary is defined by E constant, which 


along with the equation for the nozzle wall in cartesian 


2 





peusduuemtes, y' = £{x), implicitly defines the' constant a. 
Susstrtuteng for u in the velocity potential produces, 


JL 


= A cosh Nera 


from which the velocities may be determined. 





" an A гү IT, 
= = a аг —— 
E: (ir, + AS SE эх ET 





ΟΥ ОГ 
Eu ox — [m.m 
7 τι + το 2 9Υ dy 
ty 1 
а. 
p X 
9 fly + a/o)4 + x4 
por ЕСЕ DI 
эх Hy + 8/2)? * x* 
dun p Eq 2 
ay IDEN E 
0002 a NL 
3 
> = a/,)? + x2 


The constant of integration A may be determined by speci- 
fying the flow rate through the nozzle, but when the inlet 
pa Ro c i i are normalized with respect to the freestream 
ο ο ο ο is factored out of the problem. 

Velocities for compressible flow can be solved by mapping 
onthe physical coordinate system (x,y plane). 

e boundary conditions are Universal to both problems. 
These are: 


0 on the line of symmetry 


= 
u 


Ша аа at the nozzle wall 


< 
| 
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Erx) defines the nozzle boundary in terms of a ratio of the 
throat semi-height as a function of x. The throat semi-height 
is taken to be 1 for convenience. 

Pressure ratio, sound speed, and Mach number are calculated 


as before by equations 5 through 8. 
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peu DISCUSS TON OF RESULTS 


DL DNE INTEGRATION TO A STEADY STATE SOLUTION 

As stated before, the Houbolt method of integration is 
ae o all time steps. Results of the test cases bear 
this out with the larger time steps providing the most rapid 
convergence to a steady state solution. Time steps were 
tried from t = .1 tot = 100. Time steps larger than t = 100 
were not attempted because as t becomes too large the influ- 
ence that the damping and mass matrices have on the effective 
stiffness matrix becomes negligible compared to the stiffness 


matrix. That is: 


2 
SE), = S,, + 2SC,,/at“ + 6SM,,/at 


as At > œ 


Se 


SI, 


I 


The starting solutions were chosen somewhat arbitrarily. 


was chosen to satisfy the first iteration of the steady 


P3. 2At 


Solution 


S = 0 


ij?j,2at 
when the non-linear term (u) in the coefficient 


1- Mo - MÁ(y*1)u 


Beer tO zero. Ф; At and $30 were chosen as multiples of 
> > 
"uU and respectively they were 
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Ша starts procedure proved to be superior to chosing 


the first three vectors closer to the converged solution. 


J 
time steps of the previous case, the solution oscillated and 


ΤΕ 5 oat? 5 At? and oy 0 were chosen to be the last three 


converged much slower than with the starting solutions chosen 
as above. 

The stiffness, mass, and damping matrices were recalculated 
after each time step, using the under-relaxation technique 
» described above. This was necessary to utilize the special 
assembly procedures invoked by STRANL-II to prevent the in- 
admissable influence of downwind nodes from propagating up- 
stream in the supersonic region. 

Eoumocrelyecrictical tlow (M-- οσοι) and subsonic flow, 
an under-relaxation factor 9 = 1 (successive approximation) 
resulted in convergence to a Seedy Stace solution arter™only 
three time steps. Eleven time steps were required for the 
supercritical solution to converge using the same relaxation 
HI U"quersng 9 to .5 increased the rate of convergence 
and the solution achieved steady state after six time steps. 
Figure 4 compares the steady state solution for a 6$ thick 
ΙΙ. τ а ттотт ат м 2909. using the same integration 
method, with the results obtained in Ref. [5].  Chan's results 


converged in 10 iterations after using the results from the 


Bel crrtical flow as απ initial guess to the 
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ExruutIcal solution. Figure 4 is a plot of local Mach 


numbers at boundary nodes on the airfoil. 


B. CONVERGING-DIVERGING NOZZLE 
The nozzle chosen for the test cases was the two-dimen- 


Sional Oswatitsch nozzle with the boundary defined as 


A л ы: 


where the throat at x = 2.5 has a semi-height of 1. The inlet 
Mesifraken to be x - 0 and the exit was at x = 5. M = .44, 
the inlet Mach number on the nozzle center-line was chosen to 
prue ta sonic conditions in the throat. 

Two solutions were possible for this inlet condition,- the 
symmetric solution and the asymmetric solution; but neither 
solution was achieved by the finite element method. Although 
the solution converged for the subsonic case in three itera- 
tions, center-line Mach numbers deviated significantly from 
both one-dimensional theory and from Oswatitsch's approxima- 
tion [Fig. 5]. When the local Mach number M exceeded the 
ыс ас пишБет by approximately .2 (M » .64) the solution 
was invalid. Differences at the center part of the nozzle 
are due to an essentially incorrect free stream Mach number. 
Π.Ι, ΕΕ solution at x = 1.5 would improve the solution. 

A second test case was run for the supersonic section of 
the nozzle with the inlet boundary on the sonic line. The 


F oouudabvewas left tree to float. Неге the solution was 


unstable and no meaningful results were obtained. 
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IX. PROGRAM MODIFICATION 


ee element computer program for non- Lifting aır- 
us cus developed in Ref. [5], is separated into two parts. 
These have been designated STRANL-I and STRANL-II by Lockheed 
Corporation. STRANL-I generates a finite element mesh to be 
used as inputs to STRANL-II, which assembles the finite ele- 
ment equations, applies the boundary conditions, and solves 
the non-linear system of equations. Detailed descriptions 
and instructions for the use T the two programs can be found 
rer]. Only the modifications to the above programs 


will be discussed in this section. 


A. UNSTEADY EQUATIONS 
Modifications to STRANL-II to form and solve the unsteady 
finite element equations were three-fold: 


lE hcencewoelemental matrices 5С; were calculated 


and SM; 
and assembled. 


j J 


И ЗЕ бе matrices were stored on an external magnetic 
disk to be accessed and reassembled later because of 
the amount of space required to store three large 
matrices, in core memory. 

3) Me e sti ness matrix SE; апа the effective 
ια. were assembled, änd the system of 
equations solved. 

Several existing subroutines in the original STRANL-II 
program were modified to assemble the damping and mass matrices. 


These include subroutines NEWK, EMTC, DERV, and EMQT. Two new 


subroutines were added to perform the other tasks. 
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B. MODIFICATIONS TO CALCULATE THE MASS AND DAMPING MATRICES 
EMTC in the STRANL-II program calculated the elemental 


stiffness matrix by numerically integrating the equation, 


E = 7 

Equations were added to EMTC to perform the additional numer- 
ical integrations for the mass and damping matrices. All 
three matrices were calculated at the same time.  EMQT 
assembled the elemental stiffness matrices for a quadrilateral 
and trapezoidal element from the contributions of the tri- 
angular elements. Mass and damping matrices were calculated 
in the same fashion. 

Subroutine NEWK, an original subroutine in STRANL-II, 
which assembled the finite, element equations for steady flow, 
was modified to assemble E and EIE A new calling argu- 
ment, NMAT was passed to NEWK, which assembled contributions 
from the triangular, quadrilateral and trapezoidal element 


ШО ТСЕ into the global matrices S; Sep and SM; 


IS jn 


depending on NMAT being 1, 2 or 3. 
e ubr routine STORE 

Given a non-symmetric matrix stored in a banded node, 
plus the right hand side vector, subroutine STORE separates 
this system into two matrices and stores them on a magnetic 
disk. Figures 6 and 7 show the decomposition of a banded 
matrix into banded storage, and the further decomposition 
of this banded stored matrix to two smaller matrices by sub- 


routine STORE. In these figures, D, L, and U represent the 
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diagonal matrix, the lower triangular matrix, and the upper 
triangular matrix respectively. HBW is the half bandwidth 
and R is the right hand side vector. 

STORE requires an additional work area one-half the 
size of the originally dimensioned matrix which is to be 
stored. i 

2. Subroutine TIME 


Subroutine TIME integrates the system 


ο. Б πη + TUE - κ. 
by Houbolt integration. 

TIME reassembles the three matrices which were stored 
on the magnetic disk to form the effective stiffness matrix 
and the effective load vector. Once this system of equations 
is assembled, a banded equation solver is called to yield the 


solution for 4. Figure 8 is a schematic flow chart of 


J, t+HAT' 
TIME. In Fig. 8 when L=1, the lower triangular matrix and 

the diagonal of the effective stiffness matrix are formed by 
adding the appropriate contributions from the stiffness, mass 


and danping matrices. When L = 2, the upper triangular matrix 


iesteormed in like fashion. 


C. CONVERGING-DIVERGING NOZZLE 


IM Application На ОО Ето 
Regardless of the type of problem for which a set of 


system equations have been assembled, the equations will have 


Ene forn 
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READ WORK 


RE = RE 
+ UM X WORK 





SOLVE FOR 





φε-Δι 


RETURN 


Figure 8 - Flow Chart of Subroutine TIME 
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WA LEN E Is annxn matrix and Xi and К; are vectors of 
Тевес п. These equations do not take into account the known 
values of x, on the boundaries. However, for a unique solu- 
tion of the above equation, ае least one or more nodal var- 
iables must be specified and 4, must be modified to render 
it non-singular. For each equation i, either x, or R; must 
be specified but it is physically impossible to specify both 
Xi and Ri. There are a number of ways to apply the boundary 
conditions to the equations and when they are applied the 
number of equations is reduced. However, it is convenient 
ο ο μας number Of equations unchanged to avoid major 
restructuring of the computer storage. One such method is 


described below. 


If k is the subscript of the prescribed nodal variable, 


h h 


the kt? row and the kt” column of the original Kij matrix are 
set to zero, Krk is set to l and Ry is replaced by the known 
value of ху. Each of the n-1 remaining terms of R; is modi- 
r c py subtracting from it the value К: ХК. This procedure 
is repeated for all the pst values. Of course, when the 
matrix is stored in a banded mode, the algorithm will differ 
ο πα or the nmn square matrix, but the procedure is 
similar. 

Subroutine BNDRY applies the boundary conditions for 
the modified program. Setting the option parameter IOPT(4), 
in STRANL-II, equal to 1 will call BNDRY which will read the 


boundary velocities and apply them to a banded stored matrix 


by the method described above. 


σ΄ 





When the value of the χι 1s zero, as in the non- 
lifting airfoil problem, the algorithm becomes simpler than 
the above method because there is no need to either set the 


th 


k column to zero or subtract the value Ki k XL from the 


right hand vector. 
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ope Leow CASES 


ШОШ СОО Пе flow field for either the non-lifting 
airfoil or for the Laval nozzle, the following procedures 


were followed: 


1) The desired mesh was sketched with each node assigned 
a number. 


2) Appropriate input cards based on the sketch were pre- 
pared and supplied to STRANL-I to generate the data 
en punehed cards for input to STRANL-II. 

3) The above punched cards were supplied to STRANL-II 
with three additional cards as input parameters for 
each case, plus additional cards for the boundary 
velocities, if the nozzle solution is desired. 

4) Results of the finite element calculations are printed 
after each iteration, and the converged solution is 
punched on cards for possible later use. 

E T NTEGRATION TO A STEADY STATE SOLUTION 

Test cases for the integration of the unsteady transonic 
finite element equations were conducted to calculate the 
БЕА transonic flow over a 6$ thick circular arc airfoil. 
These tests were made using the same airfoil, mesh, and free- 
stream Mach numbers as Chan et al. [Ref. 5] published. These 
conditions were chosen to provide a source for comparison of 


the results. 


Freestream Mach numbers used in these calculatıons were: 


ПБД ОСЕ сеа) 
Бароу тегса) 
νυ ücnpercrugtical). 
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bach Case was treated individually with o; = 0 used as 
the initial guess for each case, whereas Chan et al.[Ref. 5] 
used zero for the initial guess for M, -.806 and then used 
the computed results as the initial guess for each subsequent 
case. 

DELT, the value for the time step, is input by a parameter 
zc scdeuxnecolumns 41-45-of the second card following the 
title card for each case when the unsteady option (IOPT(6) = 1) 
Ile Ted, 

ΓΙ Program 

a. Input 
Input cards used to generate the finite element 
mesh are listed on the next page. Cards were arranged in 
accordance with Ref. [5], in the following order: 
Mele card 
Option card 
Element cards 
Card for the total number of nodes 
Node coordinate cards 
Card for the number of boundary nodes 
Card for the boundary nodes at infinity 
Card for the nodes on the line of symmetry 
Condi for the nodes on the airfoil 
ende tor the slope of the airfoil. 
ШШЕ cards to SIRAN I cor thess tests are listed 
esse next three pages. 
b Output 
νυν... Eros TR 15 тп the form of printouts 
and punched cards. Printouts from STRANL-I are listed on the 


following eight pages. 
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INPUT TO STRANL-I 
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2. STRANL-II Program 
a. Input 
Listed on the next three pages are the input cards 
ΠΕ STANL-II program. 
ο ο τους 
Ше отр from this program is in the form ot 
Ἑ ο σουτ for each iteration and punched cards for the con- 
aoe solution. Output from STRANL-TII for M. » .909 is listed 


on the following eight pages. 


67 








INPUT TO STRANL-II 


M=0.909 


170 NODES 


— —- 


CIRCULAR АКС 


=== 


E INTEGRATION SOLUTICN 


Г СО ОА ОА СО Г СООО O +T G+ QA Tr T T O — .O —.O O in O + Gy qp Q n O м амо 
πι NNO SS ee OO πο Οσο Ωω ο ΟνΟ 
ыы ана 


CN OO XP Т O «Fr 00 47 OUYOUOIUuDnOunar-ouocouyouYOo-rgmo mro m. ОО г 
= NN MMO SF UAL DO ог Г CO 09 OS eu P OO OO eR NUI СО (СР лм ος 


— нм ам eed es = Φ oè 

O mm 

em mS 00 03 QN QS UD e P XD CONUYOUYOOOUYOf ug reoJo-9oOuoToTrTor Qe 
— NEN STL DL WD a a 00 (JM) Г оС Ооо чаа сес ОТО OO 
e = ыыы ыыы ыыы — | 

шш 

NI == 


ONDIOYDOITS WANT LIT OU On O €N VO O rO O Q г до ONON oO 
— і ο SUM OQ DONC cioe most ru 0-0 > 
-= = mm «πὶ LK 0 — — = m o ed φο 


cN Q cO JT O CAO (N QO O Q r) Q GO cO CN f CO UY OO UY COOL O^ -T O^ m cO (P çN r> — OQ CO O O 
= NNN MOM PP LO OPE DDAINDHWO OFC e eO OL O0 00 sy? «dun un Yo 
= і «πὶ e) ed ыыы 


e A A Ou en Pm 00) QS e QN + CY + OY + QO UYO QU dar G + Q Q (N ΠΟ Οι ΟΟΩώΟ ΩΩ 
NA (OA SÍ TINO OD AM CO CA GA O O —2 —=(N ANAN JY GINA O 
= — — — —— — =. — ч » о 


WONT DO T W SPO ON NPs SPP SP OO ος ΤΩΩΙ UV OU IU FUR OHO MIM MN =j = 
οι πο ώώώώωω OOP PM e O0) 0 «9 XOU Cy GA O O == cjN cc crys РЧА ОО ОО 
= 4 mi d mend md ye md mm — — — — — ου — ro! 


Е МСП СМ СО ОУ ОУ СО СМ 00400 0309 MM LADY MARINO ADONDE rco 
e сусын с MOM FFL OOP NO SPO DADDOO SSNHP De 
=i — — — — —— — — — — —— —=- ° ° 


r= r= cn Gy (n O O O YD MP (N = (rN OJCO СО rm NO TOD NE DEAN OMNMOMmO - 
NN ral el З ЗОО ОМ S OCO =A FND Ομως ών ΩΩ «ο» O 
=ч red ed med — À— — — н ed | 
Ш} 
e) 
ο ΓΟΟ O — — + Q r СОЮ (M0 Mm eO UY (Y (CO 00 (y P NOFASONOTNTONDO o 
CN m em eC rn +T + un UA O OP i Q Q O (mo cog go OO m YOU SS FAA OA OO 
= == [ne 
O Om 
O 
UY UD e СП ОУУ ССО $ O MPA MA MO O О (MDDRNLAD IS OS 090 Me СУМ e SQ e MO c 
°ч Οπή NAMATAAN S νο чс СОО O00 0 O O QOO m CIC 00 00 ST HT uu O DO 
а f een | 
11) 
LA 
WOO ND PONS MANDO AS AM ALIS ODO SM JO OUO + Y dr 00 0000 ЧМ СУМ 
CN ON o eO C0 00 Tp TF UU IO Or r> — cn P UY 4 OON QS O ο CN CON 00) 00 SF FLAN XO eC 
= — Seat St — — m ο e » 
ο NN 
O 
Q — O (N @% -+ O O +T+ QX uy (N `Q — Q — Q — — Q — +u ο ΟπιιοοΟώωώΟ -T Gy ++ QCN O O 
ем ce em — n c c ++ uy UD О ОМО со осоиш м OQOeeoJoemm TgrHruuoaor- 
© — ———— — — — — — = — 


TOP CO DUA СО О МСУ ONE © OPM p OC P HO BIA OLA DSF DNA DAMM aN 
DN =“ NAT TALO ONNO ANAD AVO NINA Mi FA A ο ο O O O 
= een ου 
un oO 
e 
OD HSHONDT TOM MANOA CJP UY — OO QY — sO AO O Т СО ГС O YN > CN = Q OUO UA © 
eN ml СОЗУ ЧАЦ МОГ СУ ЧАМ DO DO rm NIN INN) отл O 
ы” = і quand rund pud, gud gud 10 rd qmd gd rm reed peed gee Cd eed 
A UJ 
=) 


CO CO Час) ООЛО Ч COO) e LAO O) “NY O O TUNO WIN чило FDO паге QU AQ © 


UY mm NN ar ND UO Un SO SO P 09 eR CO) TP UP Q0 Un US C9 C2 QD e e CN CN 00 00) SP ST UY UY NO O O QO 


=O e ==} == een ϱ 


I 
0 


On 


68 





4 O — — CO O O O O O O O Oa SSeS eNO MN AA Se Ste © T 
ооо SGSOGS O9 999 995550555 Ο ~ 
it tk ieee Re RI  Ἱ ЕЕ — 
OOSOM-OOSOOMOSOM FSOMMOOO 9O9OnNMIFMNNMI FON MEINMNONM NO 
TOO IPFPOOmDIOlLR*ODATFTOO QQOre€uYU' doe OYOnwm TE ONODOMN un со m 
pe OO) e СС) С22 СО \СУ Сы СО ЯО ЧО © =ч сс» СО СО СОО с MMONASPODAWINAANAM AMO uu oO N 
949 ο @ @ 0 °? to 9 9 1 0? 0 9 909 99 99 99 99 9 999 9 9. 9» 96 0 0 0 ° 0 9 =— Y «ἡ 
Nt PCO MOI AN SIN SSN COP AO AO AIM (Ç (N ÇN UY e P9 (NJ 6 00) CV Uy m r (N GX rü NY O © 
ооо ооо осо ооо ооо ооо оо ооо ооо - NX OQ CQ Ο 
m ооо ооо ооо ооо оо о оо осо ооо A |x O O tN 
=i | | = 
QJ WU wW U U U U ΙΕ ΟΙ Ls LL Y LS LS US US LS LS LS UL LS LIS LL A) LI Wa Ws Lu шшш шш 
ES ONWMODOIOONYWMODIOOOODIOOWOLMNWMOOLMNIN Ouyun ONMNMONOO + ÇN 
VDOOMOOTNTOOOISOMMPOOSODONNOTEUNINNNOTMMNQANNOTUMO ο gU Ou 
ог X NNI Sud QUII. QOO T WM Ow 
| = @ ουσ 
ло mA m e e eO O00 00 Jf Un en ru e e су С С СЧ С С\З СЫ СМ СЫ СМ СМ СМ СМЧ СҮ СМ СМ ЄЧ СМ СМ) СО) πο ο 
| 
m a m OO OQO OQ OQO O m O rem meea eet Stet SS qua] (O) rl Nerd ed pod rl NS prod pu UY Ye — MT 
σοσοοσώοοοοοοοοσοσοσοοο σοοουοοοοοοοοοοοςώοσςο TF nO Ou 
| ! | ιν ΙΙ]. 111 1111. =) |= 
ο ανα Ld US UL UL LL LS LL LLS LLY LLS CLJ CAJ LLJ шыш шыш UL UL LLULLULLULLU LLULLU LU LLULLU LLULLU LLU LLI LLULLU LLL 
OOMANNA TOO OF OO DOMNNDON '-ODN<O + — r —@ c Q e eUuUMN UM un O O DN 


OO? OF ODODOND OOOW OO) 


CQ e en ON Meal M9 DN LA CA LALA Pl O O 


42 
4 
5 
0 
5 


«ео 9.9 99 9 9 99 9 999 $9 999 9 9009 99 9 99 99999 999982999 9 ml 9 ог 
Ne NIN —CNJ = — С CN, O = UN r— (N — O —4 (N (N —=ír— IE ON 0L CO CONI UY N ON EM c 
: | 

OO D= DUIOSDOIODIIDODODOIOIODDIODOIDOIIOIDOODODDIOD V UNI UN 

ооо ооо ооо оо оо ооо о ооо о о оо оо ооо ооо м СО Ф Оо 

| "ч | 
ша UO U LLL CAS LL LLL LLS UL UL LL LLS LLS LEUU UJ ULLS U LLS LL? ULI U ULLS LLS LLS LLS LLJ LLS CLS LLS LU LO (LJ ΠΕ 
Or AIM @ — O OCO O O mk ве ер «ое Гое «атаа (оТ арар Го аера Готе (ео ота» Тобо тав QN 

COCA NWOONS OD ONSWMWOD—-TONNOOMTMNNNOOSTUNNNONMTUN ο Noo- 

OY MONA DALIA PORNO OO ANAN IM SPUN OOM S COCCO с (ΠΟ Ο 

€ 9 9 9 99 9 999 90 9 9 9999 9 « e. 9? 9.9 999 * 09 09 999 9$. 99 ө 90 “л è ọm 
QD e mA СХ О ЧОЧ CN NA PNA FAAN TANN NNN NN AN NNN NNN NN NNN LO 


-- == OO "OD OOOO OOOO mira ч 
SS SS OOO OOOO OOOO O0OO 
(ү | | ! | 
LLILLJ LLULLU LLULLU LL] LLJ LL) LLI LLULLU LLILLJ LLJLLJLLJLLJ 
сз Со 2 С? © с Со Со сс © ο ο ο ο 
ошо n O O uy un O мл OO n O Сс 


А і С ні ны ыыы UCU 
ооо ооо соо O ο ο ο Όσσα Οσο 
ΓΕ 1 кг гр "H 
шшш шшш шш ш U UG US US Ld Wd DW ELI UG US U Lu LULU 
ANNO m OS JO CJ mm ooo оО eo Wm 
OMINOENOTNMLIHATOT OD OM MO A FOTON 


OMB OM BOO O =&J O UY CAO O O = GN TNN OCTAN = TOGNOLA NN O e e e O 
e * 9 ϱ 0 0? 9 9 09 ? 90 9* 9 Φ e 9 9 9 9. 99 9 Q 0 0 9 * 99 $9 9 * 99 9 9 9 ө 59 = © © 
і СУМУ (NJ CD — — NA el A — — — FOO SHO ORS A LAY NN OMA) NA A A NA MN LA ex mS OY 0g 6 [pe 


| 
QOO m -OO]QQOOQOOQOOOOQQOQOQQOOQQ]Q]QQOQOQOOQQOOQOQQ]QQOQ]QQ]QQQOQQouOogogo. MANN = 
OOOOQOQOOOOOOOOOO]QOOO]OQOOOOQOOQOQOQOOOQQOOOQOQOOOOON ος») 


1 1 


W U U ЦА CUT GL] LU U U U U U U U U U U U U U шш шш шш шшш шшш шш ш) цу 
OPNLFOONOONrSOINOOINTOOSWNOSOMWOMMOOMMONMOONOO Ў с 
eAOJU f O rN uYO Ouod Od) O ct -10) OUYOQO Qf P LA NN O O RF LAN NO OPIO ell NAPO < о 
ОО чо ОС (cum (01M OLA O «f OO eCJCOICq CJ00) SF SP UU UU XO f- P» c0 0000 ON CO 0C Q0 C0 CO CO (nO 


€ 9€ 9 0 9 8 9 5 9 99 0 9 . 9. 9 9* 9 9 e 9€ 9* 9 09 9 9 6 60 09$ 9 9 0 9$? 9 9* ὁ 9$ 9 9 9 @ @ roy 9 9€ 


AA DINA NN SY NINA FLA MOE МСО С О СМ СМ CONO NO СМ СМ СМ O (el СМ СМЧ СМ СО) an 


=| = — O m OO OOOO OOO m m rm em 
SAA OOO SO OOOO OOOO 
rie | Io) | | 
шш шшш ш шш шшш ш шшш шшш 
CQ - OO QO +T un O Oun + UG O Onn r 
сс\С^ О O OO. O + мом ооо сем 


| 
CO) O e —— — ч ч Ч tet Se IQ ч ч — .O `Q UY CO — CN O 
с) © © O OCO COCO OCO O OCO CO OCO COCO OCO OR CA OOr 
oO PI pP r T р рр í l| Y l | I 
шшш шш шш шш шш шш шш шш шш шш) LULU 
OW O I — Nui O O = OrJ =O CO O T rN Tr On O NO 
Мч оО ш СОО Ре О СОС СМ О СС СО СО ча СС) СО) СО СУ 


PID AWD aOMNOWOMD ODN POOMMANAMAAMHOMO XO CNUYm QUA QN TCO OQ P guo 
ϱ 9 9 e 9 9 9 9 09 90 9 56 O °? 9 0 060 0 0 $9 9 9 9 9 0 6 06 9 0 0 9? O 9 9? 9 °? 9 9 9 09 9€ ϱ @ = 
eA Md COO m rmm c — (NJ —— — — `Q uy FP NO ODS WIM SONA Q n ty == CN ON O US CN e m re CO) 


| 
©О=ч соо ооо ооо ооо оо ооо ооо ооо яко 
COO OOQOQ[QQOQOUQcOQOQQOQOQQQQOOOQQO«QQO(QoOOOQ]QQQOocdc«o0g]udoouocdo ч) oo 


| | e 


шшш LU WU Wd We LW LS Y Ls LJ LS U U LS LS LS LS LS Y LL UL CL] CL] CL] C.) LU LL] UU LL LU Wa Ls Ws шш 
Nor ON OFOORDMNOOMDFP OO OWININN DOWD OWMNS OMQONOO cN + 
ло NUM DALIA DNA Y Y MO O LN NN O ere YI NO O Г ~ AAN O O MA OQ ч NAO) NN PA 
VNUOWNOWMWMOMDNGOOOTIIL-OWNWO EAN NN MM TWIMNMNINOITDODOINOD © ого 
ϱ e ο 9 °? 9 Á 9 9 ο 9 9 9 9 99 e€9 9 9 99 90 9 9* 0 * 9 *.  » 99 9* 9$ 9? 9 9*9 9 0. 99 ө ° 0 
IDS wt St QICIQIN у суш СС сеу спу ШУ СУ N СЧ СЧ СУ Су СЧС СУ Су СО СЧ СЧ С eg Od eg ed Ce O Ore 


69 





e NO DMA SO 
em) FODYQM 


= = — 
С Осо ОМ 
σσ THoomM 
— = a 
DQNDNNLA Di q 
απο x — CO Ου 
=ч = met 
LA Сепо Оо 
NM A Peet DO SP 
= 9 mad 
FORD OOF Ww DO 
N QM" DO y 
= 4 ч 
MNONOOWOf 
cN CO P eU > 
= — cl 


(N 4 UON CO em UPS NU NT Oo)! 
CN OO Y UOS t= <+ ч 
— — = 


оч осмом Ф 
mat UNO e 
і ml =] 


r= NO m O mo UTA 
NN ANO Nm 
- і і 


о чо оох ОГФ 
ωσ ο LAN QA LAN 
πο і і 


u^ e m eo o C o0 «QU 
со сч сч оме 
—— —— 


TNN тг 
MOAN могил 
eod nm — = 


OO - Q" — Q ux — г 
Ou MODAN 
і чі 4 і 


more O9 
THN MOWAM 
= = ---- 


Q On O ENE rS O O T 
p — NO — — (1) Са) 4 (O 
— mi — — 


ιο οσο T üu 
Ta со Ого 
=i = = 


70 





a4 


—— ——— n -- -- 


FULL BANDWITH 





170 


CF NODES= 


D — — —À — — 


=6.CC50 


о @ @ @ G @ @ @ @ UFO 
NC. 


UP ELEMENTSs 


154 


ANC ELEMENT NCCES 


TIME INTEGRATICN TG STEACY SOLUTICN — CIRCULAR ARC ——170 NODES-—M=G.S09 


— — — — 


CONVERGENCE LIMIT 


ЕЕ. МО. 


— P аа аА 


NO. 











м — 


tt UN NN ——— MA "1 WOE (COCO f fF amr HT HT SF Hg uua uwun WOW мч ча г г =H ™ c u r Qu y tO U` UN UO а — rry piru fy Of PM atupu uu 
I ; nn ı 


l 
| | | 
ог со суг г-у одус л О {ОШ onm f 


UD CU GD UY ye 2 OS y y A A Y UD UN UP 
— — — — e ! 


NOSE 


г ло г + 


| 
Lumbre ас вым 
Р | 
съем Imo m ede «ro sr us ΓΤ изг п О асу ого сэ б {суг mo nnm. nr eo onmm.r 
Y πο ο 00 UY ty Us 


Dem nim (eo eoo оо“ mo mrorhonlowmas abdeaolrams nor 
AN AAN tet η, ΣΤ 
| | | i | 
i | 
rn 


SED [= m S meee JOU UW et A or anders oro «смс у оге ар сому усу WIND AEF TOI DIE OU WW IN 
une ANN an MT a Diu Eo FPFF т тбимгигили мила τ... a m gummmurun 
! 
Y 


26 


| 
| | 
CN mm UA go g O AO un OU VINOD AHWDIND ης DU 0) A UNO 


c 0-0 rS P ισπ σσ] UNO DO AO qao 
— — —4 —4 — — d 


> e cn | mt ad M 


+ + фллачини ол помо OÓ D DO 


an 


| 


! | 
| І 
| 
i ~ i i 
Γι | EE | 
| | | 
! ! 
| | 
) | | | 
' | | i | 
| | i : 
! | i . : | | ' 
' I | | : j А 
1 1 } | 
| ! | | ! | | | 
' | | | | ' 1 
| | | | | | i 
: ' 
! l j | ; і 1 | I 
| | | | | | 
1 ! | | { | ! | 
| | | no "ali. | 
а | | | | ' ' ' | | 
1 i | i 
| | | | | 
Em : | | | E 
| | | | | | 
| | | : | | 
| | | m | us | | 
| ! | | | | | | 
| | | | i | i l 
а | ] | i | | | | | | 
| ' | ' i | ! i 
| | i | | | { | | 
| : t | : 
o EN 
' l | | К | 
Urn ο πώ Πω Xx uo ribou unbon ο. "(0 pu YUU O учуп исэ аме учга yum Ца чо гоо О есуу у Р о ГР У С М А аг Ар Чр ОЧА очо ао фач 


Е 


l | | | 

Y || asi 

амчу чо г 9 O- nim sun ROYO rim t e пог OLD anos 
р e a et eK IN QUIN oF em «ΤΣ 


||| 


| ae i 
} 
1 | 
OM PROM ас ог о С^ О суту nomo О ус уу Ог coup Dnm rr wor- 
TT UU YU YU Q.O ο ου ο ορ κος UA UI WU 


| | | 
| 7 Poy 
| | | 


ML 





- - — nn nn, 


| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
! 


| 


uou 
οσο σῳ y 
+ — 


| 


mr 
очо сә Ох 
— — 4 —4 


US... 
ΞΡ 
wd 


NA DIA Ur 
uno uu О 


: 


| 





IIIT ыа 
——— mimo 


tO C) ct HM I cA IP m e ot ml NI 
——— — e at at tt dt —— 


{ο Φυι ος со 
QIDIVIOOIOLY 


OJIODOOVNONO + А em 


PTV VODOU U eet 
аи ааа ааа а 1 а мч а 


| 


| 





| 


і 
| 
| 
| 


1 
| 


} 


| 
| 
| 
| 
| 
| 
| 


| 
| 
I 
i 
і 
! 
} 


mou rm y dor OR ге огл © аса 
a сус PIE 09 
A ml ed ol dl rl ol ml ml ч 


INO oO PON 
NONI с 
--«------------- 


DUO чс 


—ANINN NN 
——4 —À— — I — — 


σι) οσο 
ιο ο "οσο 
[en 


ЧОР ас 


| 
| 
| 
| 


m es 


умеет? 
-π E 


| 
чуо фо 
mass 
et 


ODOT 
CN 40908 09e 0 
ee 


DINO 
NAIM MIMO 
an 


NL н 
Μον 
[a 








j 


x 
| 
| 
| 
| 
| 
| 


ee ot 
ттт unn 
td a tt md at 


FO ORD 
жолу} {чч 
444 nA 4 


| 
VOUS - mm 
MANTIS y 
——4 —- 2 — — — 


MU гусу} 
mms TTT 


Ce ee ee ee | 


i 
DO Nas 
Nenn 
444 а med 4 





ХУСТ) 


лич им ч) 


АТАНА Т АНАПА аы 





NEw NOCE 


| 

| 
TN. 
Munus O GVV 


CLE NGCE 


Zi 
| - 











1 
$ 
| 
i 
| 
i 


ooóoo 
99200 
wu Ww 
on oo 
Or ONO 
IONDMO 
озо о 
NNN H et emd 
944 9 9 


0000005000 


I 


οσο 
0000 


шишин) 
wooo 
HNO 
WO) att 


“mrad uuu 
OMR A HON а 


(oOooodooueoo 


i 


голо is Me 
о оог офф 
— — ol, cl d d pha 4 0 


m 


| 


t 


! 


{ 
| 


odoo-oobo-d--oo- 
ooooooopoeoooooo 


ου νι ТИИТ 
noonoocoop”non”ooo 
OMOOO-0MDOOMNM- OO 
LANOHOOTDNFOOMODOO 
PUMA O) ATP ποσα ο) σοι ΟΟ 
uf ean - 00^ О р 4 — ON 


° ° é ° ο °? © 0]? 9 4 9 9.9 9 ὁ 


Οσο rb O OD O O oC Ot35 


! 
| | 
{ 


MAE | о = 
0000. 00 


—— о att at 
MOODVOVOA 


ο ου 
σος O mO 
00000 yw 
ft-CjOoo ou 
порохом 
HN NT 


949 4 9 9*9 


С 


' 
1 


О ао aes 
IOHOonuhnUu 000 


'DUOVOOVDDbDOwWA ооо 


— 
о 
πω μισο Шил 
о 
e 


Y NODO 
есиме) 


€ 


πιο ουν 


“NO OVN 
силе Bar. 


rU ға) OOND Ην MIA NDW 


ОХО Dn don 


AD an 


9 9 ὁ 9 ϱ 49 9 9|9 ç é ? ϱ 9 ° 9 


OUIOOHOUHHHUOOOHONNHUIODHIDOO 


E || 
ο ὮΙ 


| | | 
| 
кыны ee w FN OoN 
AUNT F FF FANANN N 
— 4 0 «А dl A == A reo 0 


кеки 
шы 


сыс ом ао | оччо ич 
ah Ft ed hed NN NN NN 


| | 
1 | 


UY DM rN rn 
ο Νοε 
——————— 

| | 


DHOOM | 
NO my evene 


| 
| | 


| 


Л 





| 

AAA daa 
Ωω) ΟΌ 
= 


ΟΟΟΟ2ΟΟ 


0000000 


f 


O0000000PM0000000b0000000o 
Q-QOU-D0OOt20QU.foOogco 
NOITWIIOILDODOSODIMNONUN 
tC (N r Nr ry A Y UN DU PA асмо ге 
ANNDINNNIDINNNININ 
* 9484999 9291/9999 9999€ 


AHA 4 — — 4 
Odouooo 

UU UU UO UI 
OOOOOOO 


= en 
QU ug uoo 
E ' 
: | 
| 
| | 
1 >g: ον 


rame m.s FF 
| | 


E 


een 
OOOIOVOODODOOIODOODO 


μοι ο ο ο ο ο Ou 
ооооооорооооооо 
NDOMDONODR FONOY DO 
NOWNOYNPLANI+DORAO 
WONOMNONAD OMOO OMN-4 
aN A AN AAA HHA tet 
e 9*4 999 97919 9 9% ϱ Φ 9 Φ ὁ 
730000oO0O0Dboooocdcoo 


; 


ια ο ο CUA μαι 


| 


[dee 
OUOOOOuU 


+ ОООО 


оччо 4 0-4 


€ e q ..... 
0002000 


el rd ml pu rd mr 
OOUOOOU 


as Una! 


<> 
салу cO Pn enun drm em 


— m —4 — C) O) = 
QOU0occdoco — 
Wut u lU! ишшшшш ши 
0000000pP0000000 
OOODNF-OONF-MOODMO 
OOMOSOOMNDODTRN 


OCO rn: O QO 
GW OY O O (OC )C)0u) -U rs10D r. 
FAUVUNOVUO €Omv. 
со са су $ + са килил e 
Φ 0ο 0 0 ο 0 06 9ϱ|ἱ06 Φ Φ 0 ο $ о o] ç Q Q ο ο 9ο 


πο εως) y 
QT ur) ou 


€*Q00O0tO0QUnbooOdooOcdodt» 000009000 


t | 
{ | 
+104090 00 гч 
WITT 
| 
| 


] i 


curn + uy Or co 
Uu Uu uuu 


7 


ж 


| 


i 
Q^ CO Q^ O 0.0 f- O^ mno r- 70 CO 
CO C (NJ P0 PO en ta NI N d =ч 


| 
| 


ло co Q О асс 
мою омм мг 


| 


m2 DOR MIN 
NAAR OOS Eins ASS 
94 999999 


0000200 


. 
00000 
OQ OOO 


UJ ttu 
0000 
ONONO 
ONOMO 
ODIO 


ANA 00 


0000000 


| 


ο ANA aR HH NAH Het 
оч оссо) 20 U00t50u0 


JUAS US ULLI LEJ Ra) ti S U IS ES LAICA 
OOOOOOUOHODOODUIOGDOPYPOODOLUHNIDODO 
Qr-0y10]0co000 (00 00 0O0.-r 

YO UD eU) Ot) UU 

Серу’ У 


! 


ὍΣ 


| 
| 


ог οσο О ч сч 
Р Р Р-Р 09 U0 10 


E 


| I 
| 


! 
! 


@ e @ @ @ 0 ç 0|? ç @ ç 0 фр 59 0ο 
0000000>2000000000 


| 
-«οὈσσω-ῶσοω-οσω 
ONNIN FIIN 


| 
I 
] 
| 
| 
! 
| 


| | 
Е 


| 

оомооо-4оооо ч 
00000n 5 

} 1: I 
шшш тшшш, un unu 
QON009090900000 
оо О спо г} Гуго 
OONOMOHMMMON 
eo r- Damm NO 


€) O CO) en 10v em (0j (en 04 064 00 un et 

te 99 «99999 @ 6 @ 8 @ @ 

OOOOQOOOOOdooouo 
| 


1 


t 


| 


RAS Het HAMA et Ελ... 
00002000 a 

і 
О 
020000000000000900 
DOOVGD оҷ очу О 
U OI O O ASMENS ANN NIO 
UWO IOOOODATAAA TAN 
WAS FOUN Ф ЕЕ 


Тае Ете, 


ver... 
QNV0YX 


dd 


| 


i 
l 


| 


x 


| | 
| 


1 
[ x 
| 


FUN DARRO Nm y UN MIA © 
co 09 (9 90% cO Ü" O" O" u 


EMEN 


> Ын 


ΤΙ ш) иш 


9 0 Q ?@ 9 9 » Φ!|Φ Φ  Φ 5 ooo 


АТА ИИ 
t ç éd e ç 0 @ @ 


94 4 9 +оо ϱ φ 


[I| O mr Oo 








O CO OO O 
опоо Qu 


aud ao uu 
0000009 
(Qn ex drm 
vomm runs 
QN DON Fou 
ANAN AN 


Φ 49 4 9 ὁ Φ 9 


OOOo2000 
| 


| ! 


— cm TA TA —4— 
O OQCOC OOO 


Hund dei 
OOOOQOOO 
O Ouyuyuytr6 us 
чи см сим 
euer rrr 
ысу сысы сы сыс 


Oo OooOoooo 
| 
| | 
| 

WTO ar о 


CU UU quu 
mata | 


| 


| | 
Orme Or Oman 
Ata HINA 
= st tet ча 


| 


Паули ии 


| | 
| | 
| | 
ея 
0000000 а 
I 
шшш шш и) Он 
оаооооо ооо 
Πο ποσα σος 
OQQ Oo rniuinO-bmr- фл OO 


OOrw-f0o or-Do-.trowwno 
чет Cc А асо CIC un 4 


C29*2000o0O0[pooooooo 


[men 
i iiare 


ΠΕ ΕΤΕ ΥΜΕ 
0000000 
UID UIA YA UA UY 


0000000 
ως OOO OOO 


мем сзсз сзсз ср paro Pe Par 


К татат UNTEN y 
сысы су гуу) сы Ысым с ССГ) СМ) 


@ e @ Ф ө Ф 0 010 9. 4 ο Φ 9 9 φ 


OOOIVOOODOGOIOOO 


αν | 
| 


чоо оок DUA M OO: Or О 


Tara =ч а А А аф 
| | 


ΠΠ Ore I 

| 

| 

| 
за меу иом} 
NIN CIN у см со {угус 09 00 00 09 0 
| 
| = 





| | | 
| | 
1 i | 
| | | | 

"”“¿00N00-=000000-00000-00000000 OO 
O^OOOOOOOOO0O0O00000000075300000000 oo 





| Γι ! ' ' | 

LU LULU LULU UULULU [LIUU шшш шшш шш шшш шшш ш шыш ш чш 
οοσσοοσοώνοοοοσοσσσοσοσὄώὄοσσσοοοο ο 
ОИГО 7 Чи Yon Ч MOM FOWUNNINOOMDO оо 

A m FT m DOT FNN DAON iy .Q OO = сс б смог о oco 
гс] о спс 0 00 e04 A OO: NVA Aer JF um М ВАМУ Ое ил ON 
CON 00 1 P9 OS QJ mI) CO UU 0) OS m е A UN A rm O AMA eI) 0 O N 
8549 98 €? ә е 8 8 8 8 8 8:89 ὁ 89 @ @ @ @ 6/6 @ @ 6 C! @ 89 ο 


bDooonoooooOn»ooooooooooooocooooooo 


| | | 

| | 

4 а аф 4 md red ed Eel ed a rd reed FE weed nd ek reed a et ed wd md heed peed ed et eg et od et 
Q OOOOOOKQOQOOOXCO'20q]oOOOQOOqoodooodoooododoo 
j I 


UU UY w u u w U U U U U U U U U US US US US US SS US AS UU US 
QQO QOQ( OO OQOQOCO OO]Q(0]oOQOOOQOOOoo0ododoo23oo 
ооо ооо лилии ооо слили ооо 
MVA Y LY LY UY YUNCOS YOSOYS O OOO DIN) rm pP mm exsul tu tL CIC) 
00000 009 ft- ^ r— r— r— =— CO tO CO CO CU CO r) rD (D 3) Un U" U* UR CO O > 
νο NOS ON SINN ECCO 
@ е Ф @ е @ əƏ ее в @ 0 @ @ G G G 0 @ ç 0 86 0 @|6 Φ 9 9 0ο ο 5ου 


OQ OOnO0O0O]OOOQOOOQOOdUOOOOdowugoro0ooooooodoocdo 
i | 


| | 
OQ D f- Qut OU Oh- Qu fO Dt- QUyFe- wur fnm 
PS DV DOD DODO ADIT АЁ АТАТ АГА} ТУС ТҮСТҮ б СМ ССМ СМ o s 
{ 


1 ` 
4 


— — — —— —FF ш шш — — — ee eS πο 


| 
| ! 
| 

i 


| 

| 

| 
| | | 
О чече ЧОБ с О асч rus Or DO O n-rwuor-soggol 
OIL reef errera4uauunuuduuuxua-ooo-o0xg0oo09 (01 


A =— — —4 -— — 4 — ind ao o —4 —4 — —4 —4 4 — A A eA е8 с eg О et — = 


| | | 


| | | 


74 





E Ant ItLU 
€ 11 16 el 26 31 32 36 33 42 45 40 51 55 60 61 66 Τι 16 


NULES Al 
el 


45 


84 


4 


3 


2 
NOCES CN THE AIRFCIL 


NCCES CN TRE LINE 


89 93 68 104 111 118 125 132 139 146 153 159 164 168 


86 


75 


0.0 


aS 


шш 
20 
σου 
=> 
NO 
mme 

«5 


шша 


aon 


ug ee 





0.5000 


FACTOR 


RELAX. 


FACH 


-s -- 


- eee 


NES CF TTERATICNE у Тт 


——P— 


ÜELN 


CP 


Ες 


LMAC 


CLE 


VCIUM 


UCUM 


Pull 


E 


Y 
4 


NUI 


} I ! 
} 
a Ap Apr лем IN ххх“ «σαν TOT фот SIAFFT AF FF A RR NATTA 


2000 00009 OVO 0209 VOLVO YODO зо? ооооо DD VO QOOOOO Don 
и N EMILE TI aig) diy. Ia ITS v4 A 


ujJuuJ J'uusu) шшш шю ы UA UU us LU LG шш шш UU luu 
“NIDO HOODS DOA NTON SEND NN= NODT MNR 25 OD ADRIAN JAND 
ODPM NOS um O=O VETO NOS VUN AMPA NN σσ. UNION Out 
От СО) DOT DOC ἥ-πατςῳ NEON ООЛО ΟΝ” ND Ne ммт ун тосоор 


O DM lO MO DAD DO DD DD DINO e QUA ST DINA MD DN a MD O A AV O Dr DY A D νο Jp rN CG > >> > О 
5 0 0 0 ө е Ф ө . 9.9 0 06 ο ὃ ο ο э ө o €9 à €9 ο ο 9 ϱ 9 ο 9 οϱ e 9 e 9 0 ө * e 0 ? opn b ө n 0@ $9 ? 9 e ? 9$ o ο * @ 
Dune Du Nm DD Due O d e EAD Dr O NN mu eO DO WU a — CN UNT Ou mU m7 SMN DA m —— ` HHO 
каз, τι) кюк lll! ИЕ К НКЦ Л ЛК ΝΡ η | 
и ' 


i . Ü } 
OAN N wq =O) 00 NJ CNJ em SO 00 P7 L0 em X <> ποσο A OS IJ XO AJ n N rN 22 2 rry (N er" CO — MD ue IIR COM 0^9 x 0m ALT CANA 
QC OQO2U0O0O0DQGcUOOd(QcoOcoooc^ooOQo-cUGOoUOocoOOoOOOOUOOOU0oootoouoododgooowuoooooocczc 
ТИТО ОООО ОЕЕО ee a К р D л үн ү £ үлүү tq AN 
„чы шә ШШ LS LS U U AI U II U EI I LUC Gu LU II EAU U елли IE ιο.) 
me Pro o NL e Pe ALO ЧО ЧАМ "Сү = СЫ А» М ы et Л 4, ENT INT III IND UA OI DO N A ASS AN PA CU Q <> 
ANT IEN t ad PD OS hal OP чу = O NA Pr NA CO ON IA IT AL DO DON =ч лл тете JO et PAGO AO NUS me oT л 
xU ^3 VIVO pm Ie а VP IOVI AT Ον Vel SP UPS SP um em (AL Y LY O SO A MU YC OU mc) AT мч ΟΝ «ών w <= д 
Cu MO mmy VIRO IAS A Л de rd AS CQ) — CAJ CN 81 <? SP CU SO SP = SO C ЧЫ "= Сы DD ISO TO eT PR NUT ON Of ANNA 
94 4 40, 9 9Φ ο 6 9 ο 9 90 φ 9 Φ ϐ 9 ϱ 99 ϱ 9 ὁ 9 696 ο ο θ΄ Φ Φ 9 ὁ ὁ 9Φ 9:4 9 ὁ 9 ὁ ὁ ὁ 9,9 ὁ 0 9 ο ο 9 9 99 9 9 9 4 е 
M a am A t C N De Ne DELL EN SON У чуор DATEN rry — — ID UON, UP 0 0 (23 00 IT Pm es 
' I ' E ' | ! J i ! К ДЕН U ET l! Ι 


| | 
. | | М ста | 
QOOOCOOOCOOOCQOOoOO:20L72720000250oOU0o000o0(^0000)00722o000o022000200 DIDI 
νιν, A E A опр ΠΠ ΤΙ 
Us Y ud Li a US A u tu ee le JE su LS US UI dia ds LI a DS IL a tados νου νο 
Dew г IS ΩΩ πα Ρως ο DLE NA ADO IZ DO DO RN SAO — C “yz C) pus Z w C) > m OS ο” Ρον y aa O PAD 
fN CA TA a "YUY rS Q C) w мол LUN AMADO ων UA RE TNO UI A ANAL VU чле MA DUI S onu 
SSD DO ? 2 2 TY & pss pus YY Алле ТА SOUS Мм ес) TU nu ES PA LL CY: DEN 
== Je ese Хеш р сочат асо олса соха скс ee era ee el Try ee Se ee ee NA = 
a ee enter ea 


SDAA CASE A DO ADIN лол SEA AMO ABLALALI EA ALAW LAU SUN ELA SS A DSU QATA ÇNUN (AA YU Y CA IA 
! 


' ! | 
1 ' 1 [ 
i | | 
— qu — s “e w s arm «καὶ «καὶ αὶ «καὶ πὶ αὶ “πὶ καὶ αὶ κακῇ κακαὶ «καὶ шш} шш] Umi aen gummi оше] gud «καῇ umm mx S quang CUT узең gum ul a ‚тщ с=з] gummi] P aee Gud πὶ gummi) pam GE r тше σαὶ «καὶ ше] «καὶ 
SOIOO9OI9O0959DIOSOI-992>2IO09OI5999IVOIOO2OIIIPTDO9OO099999929OOO0O0D2DDO002O 
eee! Eee E V at oe ΕΕΙ ΕΥ 
LU LIU is — JU si ed his habe has ae du In e SF US Ka N J S UL Ue Lot tou: Judo uua m 
— (N 2" -— DAR PIDEN MADONNA O C5 AD DAS IDO ANA не лс сл DAN DONA Mo 
Оо ^о алг осы тыл эх? Муос FO PMO SYM FOODS NNO LOW TI NCN NI NAND T w CQ CO — TUS 
Q^QO rv i NLD Z" < n —— = n > kx 2 ο Ωω DO AO PAVO IVA VY VJ QA xQ t'N < JA C) — PT DAS DS гео лге о ма гак. 
Q — — Ç" A — — > ποσο πω 2 — Ç = OG 2 O = => O O O U* SCD — 22 — —; Z —  — — 2 NO mi (a yO me DD YD DOS e E 
9 @ e 9:9 $9 ϐ δ е © $9 01%9 υ $9 b $9 ә aa о @ bo эз hb a $99 0 0 ә е υ 0 9 е, 9 9 © © е te ee 0. ө 90 92 9 9 9 9 99 o 9 ϱ 5 


л Ссс Z 22 Q* u C" — — —% _x >> O CD > —* x — 2 J" CN 2* EQ Z y U" A CA CA VA Ü" CG* GA Os A POMADA АЛ ων ος Лл ων 
! t 
| | | | 
+ і i 
t 


! ' 


I 
| | x 


€————————————————————————————————————————— ш 
«365363090 Оо О о» ы = ος ο «ο эс ο 9 95 923 pp SS 995923 95 9 E33 9 FJD2SIS999597 C 
И UU LUCUS UE LEUR A A A E E Күк үзлз рз ΤΡΙ 
Y Wu yy ts ur LU US) UL les Glad e rto ea Ut Go D S UU Li аша шау JU IU Uy азиза a DLL 
N^ quer Ero Purus Pas) IN n AN n on ее оо 
Vors mt PM DD BP FANON MMA Nae Ç" PANS PORN МЧ = гы». ΤΈΣ UN INIRE N^. А 
PILI PE ΞΕ wee "ΠΡ нам томно N ` rJ 

δι me Abe D OS D N N et SPO SS АНЕ ЭЛ чит NT ONAN Net ON LEN 
РО a ооо ос ос A... ..o.. us. 
— — =m 203% 202" — — — —— Q — — TK DVT OPTS 337325102 YT D02 020003303077 3% Ὁ 

1 i } | i 

| | | | | 

| i | | 
ета meu o 
SA O о a O SS G — 9 C D: SS x ke im |e leas eS, 
un ara la aaa ad 
Jide dl Y Ls LA SU Al иоил a Ls UU UCD MD UU LI Шш. Шш шш. Шш о аа лы. 
ее ее рае оо ч 
DS > 
DLDORNN DO TNSR AO OMS MANAUS ών Νο ο η) ών м.о мою вчу пц с 
= p jp C CZ —= N N 1 `Ë Ç — — — WA = AY rry A = < N Ç FN =£ (Z < ur MOY Фоме т2л ION . oe BEN INN с 7 
e € 9 9:9 ϱ 9 ϱ 9 9 0 919 ə 0 ən © Fe Cee Fe 9 е əo io 0 0 0 0 0 0 0ἱ 0 ϱ 89 9 Φ à ого а о ооо оо о о а о ео ео о е 
ae ече yn = sin ADAN D ON 
pa : I | i I |! ΕΝ I | l | | poe Ge IIT 7 


| i | \ 
ve i А | $ 
ier O O SS O IS С 2528238 
== 2222 322072920" CC ο. — = = = = = 59 2253095399535 DZICISIOC IL > 2 
т Il ! нии О. ПУ Пи . d Ги na 
AP uL LU. Uu πο US ends du useful udis Ея шуи rs da AS U шаи — 
See FOS SDI FMD POST AN ONS FSC RAO SS TFA TIN DSS AION MOSM Ж 2... πο 35 
A^ am 7 m. | pe — Ορ πω Γον Ρο сле ооо хоть 
DEL > p OO — xp —* “Q 2 (> v CN cq > G 7? = “> — — —— ο ο ο G — DUM с S D O OO LODO Y was > 
ореола 
IDEM EE LUN... 2.9519 ον ο να. ο νο ο 6 6 66 Φ εἴν 4. « «* » » s ç 5» ὁ ὁ ὁ 8 ο 
= eeip gue Qmm puse m puo n uen vecti roti A jos om um oo rm mm QN Qt SN DA mt ed 


£ —.)] 
90 
Yu 


Er 


ОЕ 


i i 


| | | 
SUPT MAN ar mA Ur be A p rN = Sp rN OU MIN I nus eum m umm gis mus menu n meos mmo 
ων ο μυ αν ο ων σα σώσω ών ων] 070]u^7u0]72 0o 0( 0 VOV NO O OUN O O Q V a L O ND O O A a 2 
A A A E T ЛЛ ЛҮ т E l ob A 
иц; ss LS Ls de ων ο Εν E N LL LU LL Ds ts ser dd Us sr sr CU УШ o a ls ALS Loa e 
STAU OM sO rN Z A ODIN N Cy Sp Ç ~N р Y DA OW rn DO VA ποπ мло олсо 5 579 
т 220 2910004 — NHVSS > > (Z CA rg y < Tr 7 Up — A CO AAA "p — CO O" =p — <” — < — > SONNE OO: 
D x m UA Z) Ov rry <> — — D Ç < JA (N < =a a= NJ = Q) ο) ος) οι πιο) ΣΣ DO Ορ Sere MPMI WON Mra CS οί 
I TTS IT DIT as Pr" (Y i Ç r (` O C. C O tT i N AJ r ge ED DIDI ENDEN INT IIND ZOO DIT r: 
е е 0 9:9 е Ф b e» e ? 5/08 9 @ * 9$ @ 09 9.9 ϐ е Φ 6 ϱ 9 пә * 9 9 e$ @ 9 95.080 9 $9 9 0 е 09 9.9 0 @ * à е bo 95,9 9 @ ° 0. ө е 0 
EF TAN ft аа N 6 YJ (N (N = UA VY AJ EY utu e < CN Шу" ча Ыл СЫ АЛО ZT) C v ( > — — — UY >G — "y < — —— —— — — r (s а о Ар Ур ОСС 

ryt owed πρι οτι. mI ΠΠ πα epu eT EN N m4 


! ! : 
| | | 
t | к | 
eA CN ry <? PND DP DA NIN NA DD e Im гл шге о рсы FI < ος) et NI хо ло мм mn rr om ONS: 
κκ ο m o a — NGN NIN NN "N VN CNJ С пу У у у пе С OY T P.P p.p P F r eT PAINE WS VN EV wD 


E E ы == 


| | | 


. 


— — €—— € — --- 
"еы а == ч 
e ke ee 


| 
| 


76 





| | | | 
Ty p S r LS mes LD 5η Τσ σα ΝΕ στ 003 00 00 тү сүлгү {ч 070 =Ñ cT) Cr) 8/7 00 e 07 0n 7n e 70 en 0n en pm en en rnm en en e ^ pm en 98 0m e en e gà e en en emm uz SP 
созо ооо осоо QOOOOCOCOCC(OOC0/2]0uUO0NOGOQGOQGO^7CGUOQ^TOJOOGOOUGOCKUOOQOOOOQQOouoodougoocoocoooJ2czc: 
A еее ме герое тортом ү ү 
abs US Ls A) ed LS) aJu2uJ u dd ss YES Sd u 
eA mug PIPPI LOU Qm mu MS DDL πώ о 20 Νοκ ολο no GY mT AANTD INDAN 
ADOS DARA DOS NSPE A QA C. T CO O T uy — r= m лам etg eumgu-lguNc ο n a 2Q]Qmw' CUu Qm ewm] mus 
rm nn ом PANA DTD Emm 22.2 22! LRIENDADNMLIILIODLTNNNMNS NODNILEMNONDNM `_ 
SONNE AN OO mmN OT UNI meg qN Nem me mm A wwe ara INO OT ооо rN T A 7 — CIO US Od ute m T 
le 9 0 ө 9 o © ele © © oe Φ ὁ @ o o 0 Φ е * € HF o o 9 € e Φ Φ 9 9i o ὁ ὁ Φ Φ Φ ὁ 9149 9 9 o 9 5 0 οἱ 9 9 ç q 96 4 o o: q 9? o ө 0 у е ө өе э € a ο » 


Ox MO mu O нал 0 — CN) (N AN NIA NIN ^ч _= сы Сы Сы fN Nq 0I PP A.J CN (N (NX, — en 0m 0 0g CIO eA PNLONL C сеу т dq PO SP CO Ue A^ 
! EEE Nr кла 


| | | | 


V Amer ^ NUTAN NSD 70 00 0 — sO w N (€ - —— — = (N FN ON NN SPOT A N ANN Dd е ас АГУ N — 
Сео се) sao A a TTT ΤΡ ΤΟΣΟ 
et re ra pra pa dg qa cb p qud uM dV Pg 
КОО а ро еа оаа ла оа 
Ώ IR ND 1 ma I met NN et met ONG a) NY ON Q AG —— Q n ZG: v r— Z DAD DO агол чаму ED TI YN mt DOO a О 1 
мое чог Со) ЛР СЛ AS DIS Л SO OD ad A ID FS —— SG AQ YN: Y Y LA DA Y a ON UN UY яс UNUS dg e SO ὁ 
π ον <) — N I Aq A UN rir ША ыы ЭШ 40) ло Y na Y AQ aQ че YI VIT — Ct CN = <s — LT 1 CU гео чо iUt e GJ) ud emit) UD ee емер ае ry et CJ SA NT 
OYNS OU DEN LONDON XT ON р ге {тг cmo CY O ЫЗ а дача Уц) D UA CN t ` Y K.) Z РО ааг) а CN — G r= fN r= CN О ч a 
le e ө ө ө ө ө elo * э ө 9 à e ө e ° 0 q ° e 6 slo < o ө € ө 9 cei 0 c ө 9 ө Q ° o| o о ç 9 ө но о 4 9 o 9 9 ό:4 ө е ө ὁ 9 o co 0.00 0.0. 
ломом ЫМ σώοι (νεο ma mA NJ US MO ου ΠΠ; SN PN GON et ae I ON emt -- 


pg pg и 


| 


aR a τι ume deir ici i aai чн Pri cui iani mta ts 
III. DO MELLE OST EU COO ο ιεροί ο Ες ΕΕ Εν уу ο ες Ξ 
Pate tht itt pr ee hae Pda ran Price i ree ТОТИ te Dp ak be ге 
νιν, ud bascau didus aat Oud ees cal οὐ ου ο ο μυ το ο O A EEE а ве а № Е 
οσο A OO лину н в) а LS DISITDI SEEN ANA EST ωο απ PS Aq eme USUS m UP IO II IUS eC e TN ER INSCR Dt — 
Q^ TQ" o NINE N IT мз AD ANI DUDAR DIO гл A TIN Ur r= ЛАКОМ ОСМА ЛР N Ú TS a ODA q NOAA IL. 
= _ III SINE UN O er D U" —' N —— w —) ` e" — D as NO A A TA a N DAL NS 72 KI 2) — ον DUP PID ™ EE POs Se - 
ZA SB RODS —" — — Л Оччо DD DIM A ENON SOF TPA PDH OSD MID OPM > i; DEN SL —OquJ utMtMNSI4- 
EE ul... τ... «επι A A A A 


nur гоч Y DUNA DU UNO € 22002 DU DIAL A LAU AU A LU UV UN ET AND STE a PLAID I oP PL SI 
! | 
i | 


| | | 
| | | 


nn ορ» 


| | 
| | 


κο κο ο ροκ 


а ааа e 


i i 
=. = = 4 uud quaj dad rd cr EE comet cmt mel md a 4 um) vnd cá cud emt AD Es jaa mi eet m LI CI CI CI) — =a. _ 
PEIIISDIOOTIOTST7IOOOII 79 III IohIISOAIOSIOSISOIOIOOIIOOSIIOSOIIOTIIDTDIOSIDL 
Be a A ANA IAE р т от у я 1411. ΠΕ I | ! 
ui. —— лыш La san nl Veo m ad Vs vid nd mad LE LU und DE a us Ed a be U ως Lo nd a ES U aU SEE; Loud uideo Ge LIU LU S C Uus JE Gu Ll LU e 
CANDADO DDN —< — CO up < N — — A = DIN Y DIVA NO Z OD DA A DINO MO DD ARROZ IDAS NA DD ot COP. 
or D n co LOTION YO DAD RO O NINO A SA PATIO une mm DADO T NOUO Si DA. 
za N DIA PU UNOS a a cy S" P= OQ r— “D Q CN INE OO O rr CO | N =< 22 =p r" N a сил = Yt DI DN N IN I BD DI m. 
2204 FPS BSS BNO DOP QOD RGD ла AD PE CO ODO UA A et tt PPO VE MIN SO OV RU IN O DO =— uy r aa 82 — — — < =p — 
- € 9 9 a э 6 9:9 o @ o @ 9* $9 9 ὁ ο 9 9$ o % à 919 6 0 0 у Ф е © à 6 ὁ q © з 00; @ ο ϱ 600. 9 ο 019 ο 0? 9 o 0. 0 9:06 0 046 0 à Ф э of F © FF 9 əo . e : 
ο ο ο ὢ /ω ο νο ϱ Ὁ ὦ ων ιν ο) ων ο οπλο ο ο ο πο COS DO em ee BN eS OS s ==q == ÜA (Jx JA UA a — m ma ç" CN 2" — 
| 
1 К 
| | 
| | | i 
| | 
| Г. i | 
м CL) C) Nu) emi CI JCJ CJ mtm oma 
SVU99 70 900 YON S09 S90 2080 ооо осо оч (OO OQOQQ UQTQO OOQQO''2COCOCQOCOGOQOOOOCOCO*2QO2CO(- 
рю ри попа и l T j Pf jl b fol | I! f l Í! 1: a | 
L an AEG DJ HL I EN EI I LS UI UNI DE I „ышы ын UI UL I U 
22 7^ cene) — OR ND DAN IHN Dee 2 I LS I ?AIDI NE AMOS NDNDECCEN.INII NEN NINE ZANDD Ge" 
Nae ASIN OO NAT RNS INS Nee тог пл оков Э лоте соно то солт мА КБ ХЕС гн л тосоо нт 
ISIN 25-2554 — D — D AA OD ADO MIN GI NAT PH ро A a ru NON ^з mu ST eco ο ο τος 
NINA SD ος ποσο λα πο οσο Ñ — — DAN —.> — Z > N 2 n G > >> IST 00) HD 07 Ou С275 72Р 772545 — 2) — — p — eet ST DONS 
ео διό ο α α κ ο δὲν αὶ α вое 


56% ZDF DIT 25440 nen го LC cc IE quce e emo A DO) 
i а 
| i | I | | ' 
} s t ' 
| | | 
Í | ! | | | | ' 
Мт ES ES e ES EE eE ELOY aut Cu uuruuu t tN (NJ CN (N PI) ^N C470 0 4 49 0409 0990) Мм сыты “Ысы Сы AJFEN N UA J JkAP AJ ο ысмы" 
Sec 
' ! 


— 
30002200 OIIZS2II2IO00D0090.5 uer OC D OO OC DOO 9pgy 92so9o9929. 90900 Tao, 
логи Ра A A 


———---- 


! 
t 
jut. JU u D ш U 11-11 a tU)'U1:95-Dc us Laud Go Cr ου ле ως a a Li ha. rrt) 
ANF MANDO SDM THAN ZINN DOT NETZ x, οσο uA 2 < 2 `C a" == < ed NAT NING AOS DY ONSITE EW ᾗ 
ΞΡ. =) FD mu num NEIN DADIG De д: NAO ADVBDNEENADDBVAHAOE 5 м) а ЛАГ EL πος. 
DISINAADZAMADR TS AGA A ANA DON οσο ώσπου ο ον ЗА SAR OR om or ee COND PS NT 
SODIO Ne Ny DS OS ANNO AO r= g oc S O SNE ADA PA 
Ὁ e < < 0 @ о э Q à à 9 à O9 е 0 QQ à à 9 b 0 @ à @e е (19.06. € @ He e © © © © a. 916 е е à © €9 à 9:0 ө фо 9 9 o ту 9 9 à ὁ у 0... 010 * = Ф ο . 
а ОМГ OSE NS eR) ot NS rd SAT AN SLATES ALTE QO n p yn AA BPP NPP AANNNNNNN N mem Den 


i $ 1) I 1 | 11 t I! | { ' ! I 


— — —— — — — — — — — — — — — — — —— — — — — — — — — — — свя 2393256755$757- 
Ес CC mco mE SS AS SS Ss 
ΙΝ ΙΤ | I W | ПІ и! 
or AS ci u eh HF U U EU ds e de Lo sl II SE ELLE II AUF HUN I OHR ZU UL U una С. 
DEAD ISIDRO NAAA OS AIMAR σος Y > Su 2 IA ATADO 5 
ZIALIADZ A ARMANDO LD ANILLO DA Pt O SS > — OZ SD O DO A AN NN TOR DA ANS OS SG NS O DDN ST 
Сл LAU ISO Y DO e DANDO PAN OI NANO FINO MRI AT 
AISDISHTEISSTTAT STE NNITZT INTEL Dr Sum m DD DYCSDO SID DO NY eS SOT mee SST Tea DOS 
e 9* 6 à o è o ojo à à o э ? * à 0 à е Q € ὁ ὁ ὁ ὁ ὁ 4 Q 0 à оро à q е Q o Q е е Q €9 à 9 e * ете €9 ϐ ο ὁ Q 9 e 0 Ф е е э 9 ὁ 6:6 о а о о о а 
EN NUNC а quj. UN A х =. ~ 0 7 ч а ЛОЛ e — TN O DA AR Z — ο T" N Z A — w am qut emt nur) emi uut vun eund cod quai uud cuni — — — — — — OE ceed emg eg cmm und cmd жыз qun ND 9—2 


‹ 
t 


=m ^J 


! 
: 


= IP rw N TS TIN a SNP IA NAN A Oe IN ANIN N ποσο ν ы Са ы Сы ο πο INL OV AL P PEO EE A A 
OID FIO O 2424 GG OUO OX Q!U73 0Q]Quu 1) ] ] 2 0738:12«2C€3724]Q82€2-2€37212 QU U* 509 90 Q^* 920r 2090 C003 ω NA w: 
BEES trj ИЕ ttr ELT btt rtp paru 
Mum uU M Stud t Orca eed as) an E E CE Os y pus ts a TE RUE tque ipei рр ET 271—1). 
Selm Ste, T m CD em un N uu RU ^ MO CO) NP cem UN (NL TS e UP WA. дое хт моме ооло: ANCHE : Aqu == 
πώ πο ASE DE SL т AR as ss asas, AP ас 
Dh) ce P^ RR ОА а ла О ЛМ Реч гч м0 а О un C ui AO ата UND DN 
NUN Du = 4 Сы» wu n oun mim uu TOR V OU uut) „ымыр I Sd т т иж мм лом м ο παν жит 
€ 6€ € 8 à à э ©! 69 9 9 *9 à à à $9 9 0 0 à à» 010 ὁ ὁ ϐ ὁ ово 0 4 ϐ e е 9 09:9 9 a 9 ὁ @ 9*9 9 0 е 06 9 à ὁ b à à à à е е 9 9? 9:9 4 е ὁ 0. е e 
N SNN NT Zee ATI NZ ne N ZN ολ ων OUS Vm OO QUSE P d uA NP p σσ ο ο -- OI NOM Ad e eme Pe t — o — — LST 
Fog 1 1 | I πρ ЧИК U Es) БАШ Н О ΩΩ ΕΕ ЕИ Ир те 


| | i ' И i 

: | | 
-— 1. PUN VQ AS Z; соя; сл РГ — — оО н мер ә н тало анага ок осо асч лгал м оло NN i DA — CA = (N TADA DIAS A 
< чача ч С З Омо а о Е чч ты мы Гэ {лә чө чө CO TR CL CU WU QD CO CO dz Q AA U. ο οσο σος ος) O — —< —< — — — — — — SIN NU NIN OE OE NN TEENIE = 


а Фб - «αὶ cum) ong UM s SUM qum s j wq q j w q qum sume qum Ó s — q s s sumd cem VEM quan 4 — == 


4 


77 





TS Aq CAN Xp CO Xr UU ωσοσοσωο wn men em em en up meon en m men men m mmm 
22202002 0000009072020 20 0200 0000000000000 
ОГОО ОООО ИВТ Т О ра 
ae md DD Las LS LES U ESEL SU Y dl ci us (LJ ud LL C SL OU OL UL WU 
> OS + Den et NND EHEN O > O n An лонот лсо 
ЗУЛ ОЧ р иа t C) N — [N Us Us rn > mu in D m eY rn rN A ZO '— uy Q (Na "N 
Y NY Y Ú DCN rN I IA mm SONO OS er > `N p O Z — DINNER NAS mem OM 
DRON TENTAT CN — SD AG rN CN ry ry TQ y — QN O 22 uy A p о O yn DO) A ее 
. 54 9 9 9 9 e, € € ө ө © в Üi o о ооо фо 9 5} o ө ө 6 506 ө ο. 
I ου ο ο ο ποπ ο σπα 
ШЕ) | ‚11 ик! 


od eat mt ad NY — е а — е а “Να OL A ON O а = СЫ СЫ а — AY 
2202024920722 209092122200937223072090000232680z20u 
NOA Ee a ΤΕ er AA VI 
ed del dd id cil ed Ll 22 ш Шил UE λος, — cas Vd ead la a o мыз мы vs t а ДЈ ә а ds 
VD LD NA LO Z A Сы S2 CO OP YT AO AD e uw gd) 
ect ct SPN et Pe I re US oP POD 
Tee Pe re ЛО ee yt 
POV ME SS FT DN SP PP et LN MY DU mt 
€@ 0ο 49 4 949 4 o @ @ Q ç$ 9 e 0ο G: q Φ φ 4 9 φ °? q q ο q 9? 0@ 9 0ο ὁ 0 0 ө э э е 
"Yr (N CN emm LA ST PY ON Jmm e CP PRI NJ TN — ἀπ —s CN pa= med DP dF cet ees ee CN m UN 
uU TE o EOM He 150 S por qp { | 


\ | 
f | 
7—6 аф ted us c cmd us quuni umi qud ct ca quuni ammi a rre „ы ceed au w w umi a $ w s w s — emm 
2927272200 2000 072203Q]OCOQ237000Q]QMUUJOcJ2Q]QdQo2cocqoQqoo-JJ 
Pee Re poto er oom πι ПО τει 
~ мш. Ји Un одр. чы әм WII LL ci! dj 
De IL О СКО т у m AQ [0 — O Z rn — — r= — —— — DENT ON IMS 
TAR A NIN INE NOP MAN SD A A MA NAAA πια 
LEN TDI D UAM IO ADO NOD) A UA .5 — IMA ADO DATA SA DO r u 
rn y Y OO TEN LLC SO mS OS I SN ee Y 
«9 5 9 9 4 ' ὁ : Φ 8. 4 ϱ о фо ооо €* 9 ө $9.0 ө 06 9 ὁ Q ө ө ө ө 9 9 € + 
7 E FUR INL PP PSA PS FRITS ALA LAA AS A AO OOO DS 


. 
. 


! 
| 
! i 
| i 
) — Ç CO — I) 2 C CC) CI — uA a m mai md gnat — i — шш} ашы} cuti ci ne  — md 
SDOOIOPZSOSOIIYYo9O2I:SS7I792999090952>2y I II DIDI 9002 
ει yA ПЕЕ ОТЕЛЕ ЕЕЕ О ЕТУИ 
ICD Uu DI Iu Ioui NM da u.c оошо ү; 
οςσΦΌις м οσο σσ λσώωοσ- а а Ол о М г 
EN ed NF ODIN A PA O Y DUNA A IND ET 
«« ως AN A DNDSNWNDNN ONE WW Z5) X N N COOO Xy T = 
оно O es a DO IO =Y C: — DY DO Q — rurNÑ > — Qn e O 
9 ϱ Ὁ ө 99 e 9,9 ө 9* 08 Φ 9 6! Φ ç € ө © Q ө 0! ө 95 0 0 9 o G oe 
πετ ΟΛ ` — —— i a — N PD POR DDD DD Dee > — wd 
: { 
l i 
р ыы ed a A CS al ч (U w mah w πὶ «πὶ suns) sena CERE κκ 
30010 QOQOO0OO2OOOOOQ(qQOdqoOdcOOoQOcQoOQouoocQqoQoodocoo-oQ: 
I! ι 1 ТЕО ТЕО ЕЕ ЕЕ ОО ЕЕЕ 
RD PANS DP IA FNP OA OD ot BPM AM FAS mE ποσο ee 
ΣΙΝ юн оным млн: 
+ NDNANZMADOISAS PNDN-NDTrAUNT NT NDIDDT 
J a > 2 2 rN t'u — + ODN =D OINI NAN NCO G —— SUN AED OMT 


a $ oe 9 o o tio o q ç G э q 0) 9 Φ o 9$ € 0ο 9,9 G 0ϱ 0 0 0 © $ 9 9 9? ө ө ө 


) 
4 


е С л Сл Сао рар ао со тооло 


| | | 
| 
! i 


! | 
ey enm en 0n ч“ AQ ^" N N (N Ar ЈА) еМ М) еме) 


OSO ым м OSI ^^ еее мм 
ооо 
ООО TTE EN TEL μυ КК Л | í T үр 
US AU UL шаа U e LL A рано 
ЧО =Y T ¿p Т О ЛАС περ A FN — > — — w — > A O w Z — 2 — — < 
σολ οσοι DNS мало тэ’ ον νο σ]- ο οὓς, 
DDO DOA TPN PT Oe PANT Dr ONE ASL ENG 
— Y YO T rN AA p fN sr DANS ` AJ fN C: — — — AND EMA PO NSE > CN N — = J 
q @ 0ο 0 @ о о 0] 00 9 Φ « G οϱ ο iq Φ 0 ө @ 09 9 6:4 0 9 « Fe 9*9 9 9 06 9 o ө ө 
a a US DOT a E лс оге 2р Ὅσ рл E 
I UELUT URN NI A A II ү 
| I 
| | } 
| 
j | 
TZIDTICc ος ο Io DI IT DIT mt ә 
= = > = == =m = 


i NI 1000 
ое ны. 
ыл ы = аш DA LAR INN TNA DT DIM O — — Y DN 09-0 NS SAD a | 
E NAAA IND Zu Ze N y a O ON Ne aS TS ON ONS 
DJ a mmu INN TT IR O F NEAN FA MANA S O m a наг a 
————Q ON емен DDD cu et D ODTOD DODOL DONT ο Ὁ 2 ONA 
° € q е 9 9 4 9:9 9 9 9 8 6ο 9 9.4 9 Φ ο 9 9 о оір 9 ο 9 9 9 9 9» 94 9ο » 5 
— — —4— — — — — qui cmd quud 4 4 өф эше] ый cmi es eed ed ce ced ed oa NSN ce aud IS C8 DIN CUN `. 


i 
| 


| 
— meme mr oU 10m mem oe сым сы бы сыт Су сыы ы m AN rN (N u u ru N C u CN 
472 420020 OO n3 OO CIO DOG DOC o1 3927207 uU OO Cw 

EEE EN ET EL IR a Pt Pa 

ο το ο ο λος οπου στα ИОА ЕТЕ Ор оа рата 
AG AA VA О DDN орлы ING 
DNDN TON ώς Ο.ΟΦΉΝΟ5 ώς NAT D NN DIENT 
N DAAD ARE NENNT Du p C VA Om APD DOO AMD Ow 
2 NS me NG ОСЧ О е ао оц AS U лама Μο Σον 
. 00 0000 elo < e 9 * ϱ 9 09 9 q G q 9 0 0/06 9 0 96 9 9 9 eer woe eee 
"Φως I I NM у I gm RALIS UN Gum CDS AQ eMe ТОЛУ Сы (NJ — Г "^^ сысы = сем, 

р + | | 

1 


> < x T.O = GÇ 5 mt ÇN w — DS SSE MT NO DO DO cr NU DAA DUO 
ще лб EY > p SPP T Sp SPA p p JA VA A F.A In UA T ANAND SL u 2 2 Ç = 


— — w s s r s «αὶ «αὶ s P Lll m s s шы] ih 


78 


-- = — = ————— 


RÀ 


—————————————————— 





B. CONVERGING-DIVERGING NOZZLE 

Two test cases were run for the converging-diverging 
nozzle. The first case was for symmetric flow designed to 
ὗ E oli C conditions un the throat of the nozzle. The 
ΙΕ dealt with supersonic flow in the diverging sec- 
tion by starting with the sonic line as the boundary of the 
nozzle mesh. Oswatitisch's two-dimensional nozzle [Ref. ο]. 
y= 1 +,f/2(x-2.5)% with a semi-throat height of 1 at 
as used in both cases. Boundaries for the subsonic 
посе мете ат х = 0 апа х = 5. 

ШЕСЕ Solution 

This problem was analyzed using the mesh shown in 


Fig. 10, which consists of 126 elements and 152 nodes. In 


Bes econd card (the option card) IOPT(4) 1 indicates that 
non-zero boundary velocities at the inlet and the exit will 
be read and applied by subroutine BNDRY. This option requires 
tnat the number of inlet and exit boundary nodes be specified 
in columns 36-40 of the next card. Perturbation velocities 
at the boundary nodes follow on the subsequent four cards. 
Subroutine BNDRY reads u and v respectively for the first 
boundary node and then continues reading u and v for each 
wet aad then each exit node in the order specified on the 
appropriate card. 

S e one Case-- Diverging Section 


The mesh used for this case is sketched in Fig. 11 and 


wee ardo to STRANL=11 follow on the next page. For this 


PS 








case the options in effect are IOPT(4) = 1 and IOPT(5) - 1 
which cause non-zero boundary velocities to be applied to 


pue sonicne oniy. 


80 





INPUT TO STRANL-II 


M=.375 


NGZZEE 3 


aoe = 
оссе 
rite 
του 
момо 
πο σης) 
NO NO 
co... 9 
ANAC 
р | 

NAIA 
Sooo 
ıı ıı 
шшшш 
maths 
SOM 
eYNOI o 
99 è ¢@ 
QOO unco un 
l I! !| | 
NAN 
SOS 
| ! ! | 
шшшш 
rear 
NONO 
NANA 
99 9 9 


ONAN 


MINIMO 
O OOO 
ПБ 
шило 
O ANO 
ОЧУ СЧ 
омог 
o è è ο 
N 
ка! 
SQ CV m Nr 
πιω O OO 
ΙΙ 
шшшш 
+T 
a N 
-O mO 
O e° ° 9 e 
ΩΩ MT ON 
= | 


Y YD O DO MIA - AO O Or — nuy > A Q O r+ MNAO ONTO 
LO a s A SS SAS OSO A 
= == —— == —— —— 


cO uy r= QY O O CN +T QO UN r— CX — © С ОСО г NAMNTODDIOTEIMIN 
и 
=l = = == po) — = = =i = 


ANT SONT OVD ~AN MNT ODOT ANNT ONON SAMAN 
ο τος ο ον ο ο ONIS AOS 
— == = == = 


vx O O Ча пуке СОО AUP MAE ОЛ HONT OAND NANT ODO 
Pa Eaa T E == rb 
=; — — 4 - = —-— = 


sD O O (N ил г бл Q O CN SI COUP Qe OD CEST «QUAS QN CO CON Y «O0 09 
Ро 
— == =ч e! 4 =] =j = і 


Ιω D mt COONS ODOR DRONE OD OD SOLA POD ON SOW 
aD оет oos 
= =l = m] pd == = m -Λ 


που сс W W O U (9714 1 SPP AO OO CO CN 4 (11tí Y= ο ο ο sT' (Ty tr UY 
Hyde nz СОО ССИ Онч MOMAN = 
d = = = =l = = 


OONFAIIM NH ONLONM- N-HANMNLODPNHMNTODO 
ΗΟ DQOANNNDS — (GOO CC NGA O mi O сос сч 
— = — =- = =~ = 


-O mO m 
OQ OO 
| | | 
LLILLJLLJ LLILLJ 
лох 
чос с m 
TENN OF 
99 9 Φ9 9 
O A) < 


AAS 
ооо 
IR 
LLJLLILLJ 
ее со 
uu cO 

оом ~ 

99 9 9 9 


оо хо = 


— О O ч 
QOG(QOc2 
| | | 
UJ LLJ LLI LLJ LLI 
г 
coco cO CO - 
CN CON SO e3 Q 
оо ооо 


SM сс 


— 0 
ooo 
Fg 
шшш 
кс 
UU cO 

DORA 

eee? 


ODDO 


O O CN P" «OQ UP U^ C0 CN «TP О Q) r9 US MAT YODO MIS O DON O em = — 
TAG) mJ OO Q^ OY NU Q9) 04 4TP C49 OO Cu un O — uA CO O GQ (fN D NADO O O O 


= == = == | | == 


„=ч | |o! |] 


шшшшш 
сч СО С 


MD ACO OD OD et COLNE OD ON REMADE NO OAS COLNE QN e о 


ΗΝ QO — о ОС ао ч Оос сч о что усло LA — O QO çN <NO 


= = = = = =i m — = ad od 


| 
SN Y NA De (QU YD DON (MULA г NND ONT MAND HONT OAND =N 
OOOO m SONS Оз TS 


81 


mimi è ọ 90 @ 


AQ) a= ON ч 


—— 
ооо 
| | 
ишш 
rm 
LALA CO 
O O po Pa = 
° ° ° ° 9 


OO90- 


"y e 

> 

I | 

ш ш 

= UN 

m m 
ONO mo 
Φ 499 ο ο 


OQOOmO 


— — 
OOO 
| Í 

шшш 
= Кес 
ΥΠΟ 


Q O.O O 


OD 

ze 

πα 

e 

œx 

ш 

> 

=] 

ο. 

Ф 

2. 

=! 

O 

cx F gr ong | =l = = E = = 

ш шшш ULI 

> COU CW 

Zz P= PO MOAI PO DPM С я ССМ О СО СО О =ч СС sp OO O (N= Quy CODON ST 

O A OSI A O AA 

(>) UNA э ° @ @ -| = = — — з=) ==] к=! == 

m" CO CN CO CN 

ооо, р! 

O о «ч ADONFORODMP-T-ANINTODOTTAINMNLODON-HMNMVDONYT 

= © Sa uc ον ο омюр омот 

u. | | —- = d —= = —— 
ш ш 

O Oo yw 

== N UNO QA — cun r OO O (rN ч сс су” О оета» оч С COUPS Q^ eO CJ "QUY CA — C 

2. DODODO -Tl-Oc:rncug rN a О QQ (N.O с $ ü O pr VU C €N — Pe mago 

С) о è è = — ==} = = = a = =) = 

Ф чото ч 

ze | 

< e N ΟΟΟΩώ сха) ОСЧ ΩΩ Det ONTO ο я ССМ РО < г Oe Yun 

Сү. Q OW: AVUNA m сл ууч ушр =ч ΤΗΕ Me Om 

— — | | — = + = — = = = == 
ш u 

> хх 

о = CY G.,O tN T OCO O QA Ωμ) O O O ÇN MAM ti y> s gO O (N Pr IAN YN WO CUP WO 

< ошо СС ООП) СООСУ ШАО я ωοωώσιαοιωο- гооо мч 

LLJ Φ 9 ο ο” = = і = =] = = e j 

== CQ TO m 

YN 








O ч СО 0 00 00000 mN O mH O m O mM O m O mH O HONDO 
| | | | | | | | | — і 
Cd yd ld иин LS UL Ls 
mat PPM AMA DOND OOO DDOND OP BSP ae eet NO 
од мыло O YO DODO ο Ὀϊοσι Οώ SL QO — -O rn CY CN O —— 
-ASOLOMOMONONONONONOMOMOTOLOTANINNLUÄANNAN 
e 9 9 9 9 9$. 9 9 90 9 9. 9 ө ө 9 0 060. 0 0 9 9*9 4. 0. Φ 9 0 9 0 9 0 Φ Φ--ἠπὶ 
> а ар агр р ар ара а ао ч 


Eee 00090909000) OO OQ c0oO0O0oO0OoOocooQoooooooootdn 

οοοοοσοοοοοοοοοοοοοοοοοοοοοοοοοσοσοὄοσσσιἴς” 
— 

U W U 1) 0) UJ UY UL U US LS UA ш US LS LY US шшш ш шшш ш шшш шш ш 

СО усо сосу сеч ече (Ус О О с оао часо чохо 00 NY a Pa Pope $ $ OO 

со Оччо NT «OQ «OQ CN CN XT у сос Сосо Оч МС з е СС СҮ САС 00 e e CN 0J OO OU 

UY UYA O O O O (N (N CO dm ST UU UY NO SO PS PSUSNOS T — +T +T COCO (r) Go OO + O 

@ ° @ 999 9 9. 99 9 99 9 9 9*9 99 009 «9 0 о фо фу ө фо ουσ. 

cant ont emt ч == СМ. СМ СМ СМ СМ СМ СМ) СМ СМ СМ СМ М СМ OI ON ONIN СМ СС) СО 00 rn rn r puytin 


— — — — — — — — — — — — — — orh e pel oat ea p e pt p pel e A p pt e I Te 
@ OOO O O O cO O O OO COO COO COCO OCO O OC O OO Οο ο Ο  ΟΟοΟ ΓΟ 
ПЕ ПЕ ув = 
LO ο αι U U U U U US TU UJ UJ U US U L ο ΑΕ ο ο ο. ο ο. 1 

DS OPO DAIL OF Be SIP TOLD Bt COLA OO GY > ос о OQ fr ООО 
TOAT ~N e 00 f 00 4Q OS uy P М М о ССО асо СУСС л соо оог с 
Qs O ÇN G° Q" Q" r— DO DN DIN CO UY QO UY 7D U^ CO O Cs t— x GX CO5 (N r ON ΟΠΟΙΟ +O 
@ 9» өө 9? @ @ Φ 9* @ $9 9 0 0 0 90 9 9 9 9 0 9 9 9 9 @ 0 © © SC 9 9 SC о 
Q cO O^ CN CO CV СОЮ CN O CN Q CNN O (NJ Q CN O (NJ Q СМ СО СМ CO NN D MOA OHA Oat eet 


0ος) ο OQ C2 O <2 < C2 C2 C2 DO DIDOODOODOOOOD DODOOD OO GO e 
ο OO O O O OCO O O O O OO OOO OC O OC O 9959959 9 9 9 99 ο 
= 
CLIJLLJ CLILLJCLIJLLJLLJ 23 LJ LS L LL LLS LLS LS CLI LS ALS CAS LET CLS ALS CL! CL) шыш шшш ш ш шшш 
гоу со го о < CN CV (NC NOV CON O QO 20 00 00 00 C0 c0 CO CO (CN ONIS f P9 P9 CO 
HOO FF ODO AN SP MM DOD YI LALA a 00 00) UU 00 00 eat eh NOVO O LAN 
LAND D OO NINA E LALA LALA О Q I— — СОУСУ — SPO DMI OO Ff 
e@ ° @ ө © ө ө ө ο 9* 9 9 9 @ o 9 9 9 9 Φ 00ο ο 9 9 9 өөө ө өө 0! 


e e em em eO OQ CONO СМ IOI ONION NO ONION ON СМ СМ СО) СО 00. e слей ФТ 


ee) CD LU 
Οοοοοοοοοοοοοοοοοοοοσώοσοοσοσοσοσοσσοώσσοσσοσο ο 
Е Г Рой 

@ GY uA OO — Q Gs O @ ÇN — G° -O G° (9 C* 2 (N — O QQ ÜQ Q) QO — OX OOO +O O — 
FOL PD CD ON PENSE NS OO POA OP PPAR PASE EE 
ENTLY BN eF eTe T T eH N e TI N T TAN ONDON wO Uu 
€ 9 9 e 9 9 9 9 9 9 9 e 9 9 9 06€ 0 9 9? 060 өө о 00 0 0 9? o э ὁ 


eO P ef ef md efe ef ef cf ce Pee < С Сч 


CQ COCO O CO O O OCO O O O O O OCO ооо о ооо ооо 
ΜῊ -οοςςσοεσοσοοοσοσοοσοσοσσοοσοσοσσσσσσσσσος + 


ww w w w uw w шшш шш шш шшш ш шш шшш шшш шшш шш шшш! 
сту сс\сту со со сы см сусу сыс см см © © со со с) со сОсо сос сысы OO 
Ооо fF OQ OIN MF XT с ОООО А М М 00 (0) UY (00 = =ч СЫС ODA 
c0, 00 O O QNIN men sT xau uuodgorm-t-0og9 --—droOodgdmmoo + 
e oe 9 ο @ @ @0 9* 06 9 96 9*9 9 0 9? 9 9? Φ 0? 0 0 9 9 Φ 9 @0 00 0 0 9 O @ @ 


mt mat et met met NN NIN СМ СМ СМ СМ СМ СМ СМ ОСМ СМ СМ) СМ CJ CO CN 00 00) 00 00 00 00) ZAS LALA 


= ч + ++ ч =ч — -ἡ - o c oc + + Asco 
ee, 970729. 9 ooo o o 5 O I оо м 
| | | | | ! | | | | | | | ! 
п ишш ш ҥш Us ш ш Lu ш ш ш щш ш ш 
о сч m ο DD PR $ Tr a 92 с CGN о GN Ὁ m~“ 
О QQ (CN UA Y ”π e c NY о о пя г ГО 
DVDO-HOTODOMTMOTOMOMOMOM OTTO DIOTOAOOOMOLHN N 
e « e 9 ө ө © о ф ө ө ө ө 9* 9 oe 0 9* 9 050 9* 9 9? 9 9? 9 6 0 0 9* 9 9 9? 9 


ODDO QOuYOuOououououwuououcougouuoucodo-oot-oc 


З2Е—01 
E 
E 


2000000000 00000000000000W0o0o0o0o00o0o00o0 н 


W U U U UU U U U U UL LS US UY LL U U LU UY LY US Ls LY UY US LY UL ο ο ο. ο.) 
M0 DNNNNNNNWNO OD 00 0 1 09 LY CO (Y (Pe Pe Pe Peje e WO 


112 120 128 136 


104 


40 48 56 64 12 80 8% 


32 


сом ST DION T NP CO C092 CO «4Q Qu unco PROC) COO Un AN NO OD — O" r= Q NT ON 


ut NN NOD CO O C rÑ cu cOn ep qm uuu uM sr muU e md T cO QO cn cn O O 
€ € e 9 * 9 ө 9.9 9? 9 eee? ° @ 0 0 0 9 9 9 99 9 9 «9 * 0 0 0 9 @ 
mt et et et et AJ CUCU OI CON СМ СМ СМ СМ) СЫ СМ СУ С СУ СУ СУ Су 00 en 00 e 00 SP NP UU 


82 


Осо FO ON 
e 00 «d «0 CO QN e 
— 


MA MOU em 
SANS DQO em 
== 


+O UO CN OTO 
“(Y TUD Nr 
= 


(O QUA rut CoO 
eet USP WOM OO 
==) 


CN QOO ST O QCN O 
NN TD) — UN O 
= 


PM AQUA mm M 
a 


C2 «2 CN OO 7 ON) 
eO SPUO e DO 
= 


DO AP A OMA 
NTANDO 
E 


TOONT 
ON ST UMP © © 
= 


A DA + Po c 
NAL Do 
==“ 


O CN OTO ON 
NOM OO 
xs 


LA Ah MOON Wes 
NICYLUVO OO 


POdDUNDA SO 
NODA ODS 
= 


ANDINA cn Os 
IIND 


AND T O VNO 
сә) чоша 


CN e [P M QUIN г 
сч софа σώος 
== ή ο 9 9 
NON 








or 
Ny 
πο 


mm 


=| ==) 


ON 
NY 
= 


uc 
Nx 
-- 


ro 
NAP 
— ==“ 


mo 
nm 
—— 


CN YY 
— — 
-~ 


NA 
—— 


NG NOZZLE M=.4388 


120 
136 
152 
RCI 


19 
35 
51 
DIVE 


== == mmm 


NG 


оғо 
AMD 
d md md OE 


> 
~AD Z 
-C О 
md — — CO 


= 
ONDO 
eO MT LJ 
rd = = Ll 


ы 
Ц 4 о 4 
сот = 
ee) 


rood 
MN TO 
rj amt pe -- 


anno 

ANT 

АШ] 
k 
N 


0 
ΞΡ 


0 
ОВЕ 12 S GI LEO 


1 
0 
4 
0 
n 


Ono 


— 
— 


85 





INPUT TO STRANL-II 


UNIFORM McSH 


ys SONIC LINE INPUT, 


STEABY” SUPERSONIC FLOW 


MN ОМО ОМОШ οσο ο HONO NO MNONO NO m O 0-00 
B Donsainnamrrnmmmonni OOOOOOOOOOOOOOOcOoOOcOoO 
| r | κ | Í | | | } 
шш шшшшшш шш шш ш шшш ашшы” ш ш Ls Ww 
DO Omt OCTAN ~D LOTMAN TO 
NO 2 MMIOISDTTTANLOMTWONMN SANONAINMNMNOWOQAMNADOOUWO 
NM MIO MUT СОО СОО СЧС MOMOLFOOOTOOD-MNF-MO LT «γω 
. 9 9 9 9 9? 9 95 9 9 9? $9 9, 09 0 9 9 9 9* 9 9 9 0 
ON SATA une Mr eu eu eUuYeSIÀOm«IdH 


LA ο) Ιώ Οσο (N APUN; (N —— (N M —-MOÓÓOOOOOn»3O0oOo 
Cue ue D MONONA m ео 
} ГРЕЕТ r 
Ww U U U L U LLS Uu UL L? U UU US uu UU UU UU LU UO LL 
+ Q CN O O Y DOPNILAS0YNDYN SO 
van, NVYWM-NONMWMOO MALA MO) > J.O О — O O Q UA OY €N CY QO O 
оо (NUYU 41 COOCOAOCNUYA JS: MINO QA 0 e 0D XT CJ4O UY OD NA LALA 
ο ο . 9* 9 9* @ 9 o 9 0 9? 9 0 0 9. € 9? 90 909 9 9 9 


OO LAI) OUND QO D +à Ca а а аа сусу Сч 


ry о ky o Soo o 
m NASAN TEMO MVU NANN roD O O O O O O OCO O O OO OCO O CO O O OQ OO 
І ДЕ ЕР | | І | 
шшш шшш ы”шш ш UU U UU U UU шыш ш Ls LW La 
FUN NO PO NO Mp DM DRNA TY MO $ De DN 
SPN WD O RMN OD HAO ώς Ос Отсос СКС ООО осме чос 
ло мал О МОМЕ М ОО Ç GN S O O (N + + O D O OQ += O (rN G CN 
° € ϱ 9 ο 9 9 9 9 °? 9 9 0 өө оо о 9 o € 
ug (YA NO NOD C0 O5 00 QN (0) e 00 0m 00 e ST e ч 


mU O N MUT CU O em) Mr «Or Un eO m mM — — — — — C) Ç) CO O O CO CO COCO COQOQ 
“Sr  ὉώΟο ορ Ωω οοςυοσοοΟοΟΟΟΟΟΟΟΟΟΟΟΟΟςΟ«Ο) 
р} πι tot 
zs MR MH C SUC EUER ERE C 
DOOM ANMOIAAMODOA FMNO OO 
со ОсОО CN d un 00 OCIO eU e CDU C9 UY NN a. an OOO mo 
OO Net MONONA ~st Ψώωωώ Ωω ώσσωροΟοώσώ ο οσοσοοώώ ο 
° 0 e 9 0. @ ө 9 9 9 9 9 0 0 0 0 9 Q o 9? 9 9 9 
OO DAMA DÍ Di MM mM mM M IM AO VOS ON 


OC NeT осо чеч го CO OO O COCO O CO CO COCO CO O CO OCO COCO Nr 
| | ııı ss ıı ı dO d 4d d d jp d gd ! | 
Ши] ujuJuluJuJjJuuJjululuJjttiujulujujtulujujul 
Ў COLD SO Mh + + O СО г < АГ ODO LO e 
OM mF CO UN Xd ur o9 O —= c.f O O CY += GOO UA — DO LA CP OO ФО —-—0) mm -Uun ho 
ET PME OH ON OND ч MIETEN MIDOOM LORD LOS TO WO 
. 9 949 9 @ @0 9 9 9 о 9 9. 9? 9, ? 9 °? 9 46 9 0 
DMN М ed Pee e I 4С яс) чс су суя 
° | 
— Q* Q^ eC «df Of бло ССОО © O m nn m0 OO mem HAI O OO DO OS SO TI O O OO «ος. 
+ ОСЧ AP rO m AO OC OCO ОО ОО О О ООО ОО ООО du 
1 1 | 
ο ο ο πο Wai il 
CON SO DOP LAA SO ü — (n DOD Ca O pun OO 
DON ML rm BO © MIAO FOR DMA SON FADO WRI хо + οσο 
оо =? οσο οοώμήώ соге о Сосо ге су ОСО СОС 4O FO ΟΟΟΟΟΩΙΩ THE 


ооо mO MN QO IA C — —— —————= NON NON 

° | 

а ССО ОСМОН с О чо ООС амф | >) ч ч о чо >) - - + <q mo 
O = pr cn O CN OC UY MM O ο O O O O O O O O MY 
| D» | | ! | | r | | | | 
ш ш ш ш ҥш ҥш ш ҥш и> и ш 
со 0 Ὁ FAN Fr + r D ч O 
DOODAMSFOMHONMNMDOAM $ A N M =— ο ο O Mun 
Quo а A oo Nm 
оо @ ..... ..”-.eaeae.„..,:— CC, €? ?* 9? 9$ 00 

πιο Y OUNNONOSNOOODODVOMOM-ODOXD 

о 


9 СЭ = cO uy O OD —4 CN) xT UY F< 5 O (NJ (YU — 3 — — — OD OOOO OO OOOO O NR 
QG) CÇGCONOCNUA AURA SIOUSTWNMYOO COCOCOC CO COCO CO COC OO OO OO OD τν 
ο μα αν ασ... 
LLJLLILLILLILLILLJ LLILLILLJLLJ LLJLLJLLJLLJLLILLILLJLLJ LULLJ 
WO FOP PMSA O ADIN FeO WIA OO 


ο MAT OM DONADO чп) TF WANN uY OOO гг A DO DIS TN O O UN O 
OO, CQgojucJuveT — ГУ пум CJUO) P Jue C10 m LM OT IL = 
° 9 99 @ 9 9 @ 0 9 9* 9 9. @ 9 9 9 9 9 9 ὁ 
оо Am UO CN] xQ UY Ç) CD — — — —— —— — N (N (N CN 


84 





ONETO 
сто Ф 


LO e MD 
mt OAT WO fh 


TOON 
— OT Uf 


е pm 
A CN T Or 


MN Ow Ow 
et QT SOT 


чес Ov us 
і МУ УГУ 


O WN U ST 
= NTA 


1.31 76-02 1,289E-01 T.T186bE-01 2:215E-O01 225T6E-OT 2v8(3ESQOT 39 T0500 


me MU 
MIM) MATAS 
. 9 
φ- (τὴ 


85 





SH 


TTT 
f/f ffi 


86 


Figure 10 - Finite Element Mesh for the Converging-Diverging Nozzle. 
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